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Abstract 



Here we correct, extend, and clarify results concerning the spin Hamiltonian 
Tis used to describe the ground manifold of Hubbard models for magnetic in- 
sulators in the presence of spin-orbit interactions. Most of our explicit results 
are for a tetragonal lattice as applied to some of the copper oxide lamellar 
systems and are obtained within the approximation that Tis consists of a sum 
of nearest-neighbor bond Hamiltonians. We consider both a "generic" model 
in which hopping takes place from one copper ion to another and a "real" 
model in which holes can hop from a copper ion to an intervening oxygen 2p 
band. Both models include orbitally-dependent direct and exchange Coulomb 
interactions involving two orbitals. Our analytic results have been confirmed 
by numerical diagonalizations for two holes occupying any of the 3d states 
and, if applicable, the oxygen 2p states. An extension of the perturbative 
scheme used by Moriya is used to obtain analytic results for TCs up to order 
(t is the matrix of hopping coefficients) for arbitrary crystal symmetry for 
both the "generic" and "real" models. With only direct orbitally-independent 
Coulomb interactions, our results reduce to Moriya's apart from some minor 
modifications. For the tetragonal case, we show to all orders in t and A, the 
spin-orbit coupling constant, that TCs isotropic in the absence of Coulomb 
exchange terms and assuming only nearest-neighbor hopping. In the pres- 
ence of Coulomb exchange, scaled by K, the anisotropy in Tis is biaxial and 
is shown to be of order Kt^X^. Even when K = 0, for systems of sufficiently 
low symmetry, the anisotropy in Tis is proportional to t^X^ when the direct 
on-site Coulomb interaction U is independent of the orbitals involved and of 
order t^X'^ otherwise. These latter results apply to the orthorhombic phase of 
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La2Cu04. 
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I. INTRODUCTION 



A long standing problem which has attracted much interest recently concerns the mecha- 
nism whereby spin-orbit interactions give rise to magnetic anisotropy in magnetic insulators. 
This subject, which was extensively investigated three decades ago,[ 0, @] has recently been 
the object of renewed attention due to interest in the lamellar copper oxide systems. [ The 
first of these to be extensively investigated, La2Cu04, has a small orthorhombic distortion 
away from a tetragonal structure [ ^ and the above mechanism was shown [ ^, H, to give 
rise to an anisotropic exchange, including that of the antisymmetric Dzyaloshinskii-Moriya 
type. In that system there are two anisotropy energies. [ H, |^ One of these, the out-of-plane 
anisotropy, is of the form aN"^, where is the z-component of the staggered magnetization, 
the 2;-axis is taken to be perpendicular to the copper oxide plane, and a is an anisotropy 
constant. This energy causes the spins to lie in the basal plane. There is also an in-plane 
anisotropy energy which selects the orientation of the spins within the basal plane. Until 
recently the discussions of the origins of anisotropy were confined to the orthorhombic struc- 
ture. However, more recently a family of copper oxide materials of similar structure, but 
which are actually tetragonal, have been studied [ 0, and found to have roughly the 



same out-of-plane anisotropy as La2Cu04. The earlier studies [ |], ||, 0] did not predict any 
anisotropy in the tetragonal limit. Accordingly, a reanalysis of anisotropy for the tetragonal 
systems ought to show a common origin of the out-of-plane anisotropy which does not rely 
on the orthorhombic distortion. That is the main purpose of this paper. However, in the 
course of this work, we have found that a number of general questions concerning both the 
results and the methodology required some clarification, which this paper is intended to 
provide. 

A microscopic basis for superexchange between magnetic ions was first given almost 



forty years ago by Anderson. [ ^ In the language of a Hubbard model, [ |13[ his calculation 
started from an orbitally nondegenerate band in which there is one electron per site in the 
limit of large Coulomb interaction U whenever two electrons occupy the same site. If the 
kinetic energy is completely neglected, each electron (or hole) may be characterized by its 
spin. When kinetic energy (described by hopping) is included perturbatively, one finds a 
spin Hamiltonian, which in low order perturbation theory can be expressed as the sum of 
contributions ?i(i, j) from each bond (i, j). This spin Hamiltonian describes the perturbative 
removal of degeneracy. In higher order in the hopping, one encounters contributions to the 
spin Hamiltonian from plaquettes (at order t'^/U^, where t is a hopping matrix element) and 
eventually from even higher order clusters. Ignoring higher-order contributions, Anderson 
obtained an isotropic exchange interaction between nearest-neighbor spins, 

n{z,j) = j{t,mz)-s{j) , (1) 

where J{i,j) = 4:tfj/U and tij is the hopping matrix element between sites i and j. 

Soon afterwards Moriya [ used Anderson's formalism to study the effect of spin-orbit 
interactions on superexchange between magnetic ions. He showed that for sufficiently low 
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symmetry the most general effective spin Hamiltonian for two spin 1/2 magnetic ions, sucfi 
as Cu"'""'", is of tlie form 

H{i,3) = J{t,mz) ■ S(j) + D(z, j) ■ S(z) X S(j) + ■ M(z, j) • S(j) , (2) 

where M(2,j) is a symmetric 3x3 tensor. The first term represents the isotropic sym- 
metric exchange. The second and third terms represent the antisymmetric and symmetric 
anisotropics, respectively. Moriya's results were obtained to second order in the hopping per- 
turbation, but in principle provided a framework in which the spin-orbit interaction could 
be included to arbitrary order. Convenient explicit results were given to lowest nontrivial 
order in the spin-orbit coupling constant, A. 

Much more recently, Thio et al. [ ^ found that La2Cu04 is described by Eq. (2). 
Consequently Coffey and co-workers, [ ^ invoked this Hamiltonian to describe the CuO 
planes in the cuprates. They found that D{i,j) can not be the same for all bonds < ij >, as 
was assumed by a number of previous authors. The form of the is determined by the 

symmetry properties of the crystal structure. The first attempt at a microscopic calculation 
of the vectors D{i,j) was made by Coffey, Rice, and Zhang [ |^ in the framework of the 
Moriya theory of the anisotropic superexchange interactions. Within this theory, D(z, j) is 
of order A, whereas M(i,j) is of order A^. Therefore, many authors neglected M. Naively, 
one expected a gap in the spin-wave spectrum due to anisotropy, and this is what one finds 
when M(2, j) is neglected. Subsequently, Shekhtman et al. (SEA) [ have shown that 
M(z,j) can never he neglected. Most interestingly, when M(i,j) is included, they found a 
hidden symmetry in 7i(i result of which, inclusion of spin-orbit interactions did not 

reduce the degeneracy of the ground state of the pair of spins [ Q Their result was 

that 7i{i,j) could be written in the following form: 

mj) = - • S(j) + D(z, j) ■ Si^) X S(j) + S(z) ■ ■ S(j) , (3) 

where the vector D{i,j) is bond-dependent and [A ® B]^,^ = A^B^. As SEA show, the 
result d^) indicates that although the pair interaction is not of the isotropic form of Eq. 
(0), it is rotationally invariant and hence the energy level spectrum of the pair interactions 
consists of a singlet and a triplet, just as it would in the absence of spin-orbit interactions. 
In previous work the terms in D(z,j) and those in M^^,^(z,j) = Df^{i, j)D,^{i, j)/{2J) were 
not treated on an equal footing, and therefore this hidden symmetry was never noticed. 
Furthermore, SEA showed that even though each individual bond might have this hidden 
symmetry, the crystal as a whole could have anisotropy because of the frustration caused 



by the competition between exchange interactions of different bonds. [ [T^ In particular, for 
La2Cu04 they found that the anisotropy was a result of this frustration. 

All the work cited so far relied on the idea, introduced by Moriya, that the effect of spin- 
orbit interactions could be taken into account by a gauge transformation on the hopping 
between sites. As used by Moriya to obtain results up to order t'^/U, this formulation is 
correct and convenient. However, this formulation does not form a correct basis for calcula- 
tions to higher order in t/U. Thus, as we shall see, the hidden symmetry of SEA, although 
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maintained at order for constant U, is broken at order for constant U or at order 
for nonconstant U. (Here constant U means that the Coulomb interaction between holes in 
two orbitals does not depend on which orbitals are involved.) In addition, the calculations 
of Shekhtman et al. [ |T^ or Bonesteel [ |rB| for the anisotropy of the cuprates were based on 
terms requiring the existence of a distortion from tetragonal symmetry. However, the easy 
plane anisotropy is observed [ |10|, |l^] to have similar magnitudes in both the orthorhombic 
and tetragonal cuprates isostructural to La2Cu04. The main reason for the failure of the 
previous calculations to give anisotropy for the tetragonal cuprates was the fact that these 
calculations neglected the Coulomb exchange interaction. From the results of Barriquand 
and Sawatzky (BS)[|T^ one can see that they partially included such interactions. However, 
it remained unclear which aspects of the BS results would persist when the calculation was 
pursued more systematically. In fact, in Ref. |T8]it was shown that for tetragonal symmetry 
Coulomb exchange interactions played a crucial role in determining the anisotropy. 

In view of the above history, the following points remained to be clarified and are ad- 
dressed in the present paper. 1) One should generalize Moriya's results for 7i(z,j) to the 
case of nonconstant U . Having done that, we find that when reduced to the case of constant 
f/, our present results differ in a small way from those of Moriya, who overlooked some 
A-dependent contributions to the magnitude of the isotropic exchange interaction, J. We 
also give general results for superexchange interactions, i. e. for the case when the copper 
ions are separated by an intervening oxygen ion. However, the results are given in a general 
form which can equally apply to systems of ions other than Cu, as long as their ground state 
is orbitally nondegenerate. 2) Since earlier calculations for the cuprates omitted hopping 
between excited states of the Cu ions, we have reanalyzed the role of symmetry at arbitrary 
order in the matrix elements tcu-Cu = ^ which describe the effective hopping between copper 
ions. We find that in the absence of Coulomb exchange interactions [i. e. for K = in Eq. 
(I), below], one recovers isotropic exchange for a simplified "generic" model which describes 
the complete 3d band for copper ions on a simple tetragonal Bravais lattice. This isotropy 
is the result of the high symmetry of the crystal field levels and the resulting high symmetry 
of the hopping matrix elements. This result shows that in the absence of Coulomb exchange 
terms, one retains isotropy in j) to all orders in both A and tij/U and thus that inclusion 
of Coulomb exchange interactions is essential to obtain anisotropic exchange interactions in 
the tetragonal case. 3) For a tetragonal lattice we find that this accidental isotropy in H{i,j) 
is removed at order t^X^K when Coulomb exchange interactions are allowed and we give 
detailed expressions for the exchange anisotropy in terms of the hopping matrix elements 
and the matrix elements of the Coulomb interaction. 4) For a crystal with arbitrarily low 
symmetry (i. e. when the crystal field states have no special symmetry), we expect to (and 
do) find a removal of degeneracy of the spin triplet. This breaking of rotational invariance 
occurs at order for the case of constant U and at order when U is nonconstant. These 
results modify the conclusion given in Ref. |^ 5) In contrast to all previous work, we also 
found an in-plane anisotropy originating from the anisotropy of the spin-wave zero-point 
energy. 
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Most of the above results have been obtained analytically, both for the "generic" model 
(with only Cu ions) and for the "real" model (in which the Cu ions are separated by oxygen 
ions). Furthermore, we have corroborated our results by comparing them to results obtained 
by numerically diagonalizing the Hamiltonian which describes all possible states of two holes 
on one bond. For the "generic" model, there are 20 single-particle orbitals, 10 on each 
copper ion, so that in all there are 190 two-hole states. For the "real" model there are 
six additional 2p states on the oxygen ion, so there are 325 two-holes states in all. In the 
tetragonal case, where we know that the exchange interaction matrix J^uihj) [see Eq. (|13[), 
below] is diagonal, its values may be deduced from the values of the energy splittings of the 
ground manifold, as is discussed in Appendix A. 

Briefly, this paper is organized as follows. In Sec. II, we first introduce "generic" and 
"real" Hamiltonians and then discuss the perturbative framework we use to calculate the 
spin Hamiltonian Tis- The actual perturbative calculations oiTis are described in Sec. Ill, 
although many of the details are relegated to Appendices. Here we give expressions for 
J{i,j), D(z,j), and M(z,j) up to order t'^X^AHc for the "generic" model and analogously 
for the "real" model, where ATic represents Coulomb interactions beyond the approximation 
in which the Coulomb exchange K is zero and U is constant. The case of tetragonal site 
symmetry is discussed in Sec. IV for both the "generic" and "real" models. There we 
prove a theorem, valid to all orders in the hopping matrix elements and spin-orbit coupling, 
which says that for nearest-neighbor hopping, the complete spin Hamiltonian is isotropic 
when Coulomb exchange is absent. Here we display explicitly the leading contribution to 
the anisotropic exchange when Coulomb exchange is treated perturbatively. In Sec. V 
we discuss the experimental consequences of these results. In particular we estimate the 
anisotropics and spin-wave gaps which our work would predict. In Sec. VI we study the 
case of arbitrarily low symmetry for the "generic" model and show that the anisotropy in 
J^u{i,j) is of order t^A^ for the case of constant U and of order t^A^ when U is not constant. 
Finally in Sec. VII we summarize the conclusions of this work. A brief summary of our 
major conclusions has been given previously. [ O, 



II. HUBBARD HAMILTONIAN AND SYMMETRY OF EXCHANGE 

A. Generic Model 

In this section we introduce a general Hamiltonian, versions of which will be studied in 



this paper. We start from the following generic model, [120] which captures the symmetries 



of the cuprates. For holes which reside only on the Cu ions, this model is given by 

i,a,a holes, h 
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Here d]^^ creates a hole in the ath spatial orbital, whose single-particle energy is eia, with 
z-component of spin a on the Cu ion at site i. In general we allow hopping with matrix 
elements tiajjj between the a orbital on site i and the (3 orbital on site j. This Hamiltonian 
also includes direct Coulomb interactions between electrons on the same site (scaled by U) 
and on different sites (scaled by V) and exchange Coulomb interactions between electrons 
on the same site (scaled by K) and on different sites (scaled by N). Our numerical work 
indicates that when t 7^ 0, the effects of V and N are not qualitatively different from those 
of U and K, respectively. Since the latter are often dominant, we shall neglect V and N. 
In principle one should also include Coulomb terms with four states, U ap-ysdl^^dlpgdi^gdiSa . 
Here we follow most of the literature and start with the simpler Eq. (^, which involves 
only the Hartree-like terms U and the simple Coulomb exchange terms K. Equation (^ 
can easily be extended to the "real" model in which we include p states on the oxygen ions, 
with hopping between them and the d states on the nearest-neighboring Cu ions. 



B. Single Site Hamiltonian for Non— Interacting Holes 

In this subsection we briefly discuss the basis states used in the perturbative scheme 
described in the next subsection. We first consider ions within a single-particle picture. 
We therefore start by considering the effects of the crystal field Hamiltonian, Tij., and the 
spin-orbit interaction, Tiso- The former is constructed so as to give the observed ionic levels. 
Including only such energies the single-particle Hamiltonian is 

/3))^^,rfLrf,^., , (5) 

iaa ia/3 

where 

era I ^ era 

in which < ia \ i(5 >= L^^ is the matrix element of the yU-component of the orbital 
angular momentum between the two single-particle states, and is the Pauli matrix. We 
shall often present results for ^-independent matrix elements of L. 
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For many purposes it is convenient to diagonalize the single-particle, single-site Hamilto- 
nian Tix + Ti-so- We may choose the wave functions \ ia > to be real, in which case the matrix 
elements L(^^ are purely imaginary. As a result, every single-particle energy of Hx + Tiso is 
at least doubly degenerate. That is, the two linearly-independent wave functions which are 
related to one another by time reversal, 

V'a = "^{Vaa \ « T> +Zaa \ « i>), 
a 

0a = E(-<a I « T> +yaa I « i>), (7) 
a 

belong to the same energy. We use Greek indices to label the crystal field states in the 
absence of spin-orbit interactions and Roman ones for the eigenstates of Ti^ + 'Hso- The 
latter can be characterized by pseudospin quantum numbers, a = ±1, and are associated 
with the creation operators c\^„. These operators are related to the d\ 



^iaa ~ $Z(^Q!Ci)o'o-i'^m(Tii (8) 

where the unitary matrix m^^ is 

<a = i ~ [f ".f ) = <al + ^<a " ^ , (9) 

where is a real scalar, v^^ is a real vector and I is the 2x2 unit matrix. This leads to 
where 

(Here the dagger operation on m^^ operates only in terms of the two by two matrices 
as in Eq. and is not to be applied to the scripts z, a, or a.) The transformation 
that diagonalizes the single-particle Hamiltonian is in general different for each site. 
Consequently, the single-site energies may depend on the site index. However, in certain 
situations, for example, in the presence of the tetragonal to orthorhombic distortion in 
La2Cu04, it is possible to define the transformation such that the single-particle energies 
are site independent. This will be the case for some of the explicit calculations which are 
presented below for the cuprates. 



C. Formulation of Perturbation Theory 

For Cu^"*" ions in a configuration we are dealing with an ionic ground state having 
one 3(i hole whose spin is arbitrary. When we include the oxygen ions in the model, those 
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ions have filled 2p bands in their ground state. In either case, in the absence of hopping, i. 
e. for t = 0, the many-electron ground state manifold is one in which one hole of arbitrary 
spin resides on each copper ion. The energy levels within this ground manifold, when the 
remaining terms in the Hamiltonian, especially hopping, are considered, are the object of 
our study. 

When hopping is introduced as a perturbation, the splitting of the hitherto degenerate 
ground state manifold can be described by a spin Hamiltonian, Tis- In view of time reversal 
invariance Tis will consist of two-spin interactions (between nearest and further neighbors), 
four-spin interactions, and so forth. In the present paper most of our results will be for the 
nearest-neighbor two-spin coupling constants, except for the general theorem of Sec. IV, 
which makes no assumptions about the specific form of Hs- If we only consider two spin 
interactions between nearest-neighboring spins, we effectively write 

ns = Y.nhj), (12) 

where {ij) indicates a sum over pairs of nearest-neighboring sites and for spins 1/2 

n{t,j) = Y.j,,it,j)s,it)s,{j) , (13) 

where /i and z/ label Cartesian components. We refer to the case when J^^uihj) = J{hj)^fi,v, 
where 6 is the Kronecker delta function, as isotropic exchange. [To avoid confusion between 
the two kinds of exchange, the terms in Eq. proportional to K are referred to as Coulomb 
exchange.] Appendix A contains a discussion of the possible anisotropics in TC{i,j). 

The major objective of this paper is to discuss the symmetry of the matrix J{i,j) and 
develop perturbative expressions for it on the basis of the generic Hamiltonian of Eq. 
and its generalization to include the intervening oxygen ions. /^From our point of view the 
most important early work was that of Moriya,[ 0] who studied a simplified version of the 
above model. The most significant simplifications necessary to obtain Moriya's main result 
were to neglect the Coulomb exchange, K, and to assume constant U, i. e. to assume that 
Uia^ip did not depend on either the site index i or the orbital indices a and /?. In particular, 
when Uia,ij3 is independent of a and /3, the wave functions for the two-hole states are Slater 
determinants of the one-hole states as obtained by the canonical transformation of Eq. (^. 
In other words, in this very special case, the exact eigenstates of the Hamiltonian Ti^ + ?iso 
also diagonalize the Coulomb interaction. Tic- In terms of these new single-particle states 
the transformed hopping Hamiltonian now assumes the form 

■id 

where 

^^. = E01),,,4.c,w (14b) 

ab 
crcr' 
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represents hops from site j to site i, and (t^^^,) is the 2x2 matrix 

(tl) ^ E i-,./3«a) = + ■ a , (14c) 

a/3 

in which A'^f^ (^ab) ^ ^^^^ scalar (vector), that can be found using Eq. (|^) and the 
representation in which tiajp is reaL By hermiticity these coefficients obey 

Actually (and this seems to have caused much subsequent confusion), Moriya did not write 
down Eqs. ( |14a| )-( [n^ ). Instead, in a further simplification, he truncated t to include 
only hopping between the t = ground states. Even for his calculations at order t^, this 
simplification is slightly incorrect. However, we should emphasize that this truncation is 
totally inappropriate for a discussion of effects of order higher than t^, since hopping between 
excited states then comes into play. Also, when U is not constant, hopping between exact 
eigenstates of Hx + Tiso + 'He is no longer a single-particle interaction. To see this note that 
there are matrix elements between an initial state, in which both holes are in their ground 
states on different ions, and a final state in which, for instance, both holes are in excited 
states of one ion. Such a process explicitly relies on the fact that the two-hole states are not 
simply obtained from single-hole states. Thus, in this case, when "final-state interactions" 
are present, the hopping perturbation involves four electron operators. 

Accordingly, to study the case when U is not constant and when Coulomb exchange is 
not neglected, we write He = Heo + A7-^c, where 

'HcO = -^Uq'^ dleeadlfSa'd-i/Sa'diaa = ^ ^0 '^laa'^lha' '^iba' Cjag- , (16) 

icK/3 iab 

ctct' era' 

and the additional Coulomb terms resulting from nonconstant U and K take the form 

i aba'b' 

with 

ace' 

Kaa'a^a',{^abb'a')=J2Ka'(i^aay^a'a') , ((mU)Vl,,) , . (18b) 

aa' 

Expressions for AU and K for tetragonal crystal-field states in terms of Racah parameters 
are given in Appendix B. 

In the following, we will calculate the effective spin Hamiltonian using perturbation 
theory in which we take the unperturbed Hamiltonian to be 
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'Ho — Tix + '^^so + 'H-cO — X! E^iacl^^Ciaa + q ^ clao-cJ^o-'Cifto-'QaCT (19) 

iaa iab 
aa' 

and the perturbation to be 

V = T^hop + ^n, , (20) 

where these quantities are given in Eqs. ([1^) and (|1^. 



III. PERTURBATIVE CONTRIBUTIONS TO W(I,J) 

A. Contributions of Order t^ 

The lowest order contributions to H(i, j) are second order in t. At this order in t in the 
absence of the Coulombic perturbation A?ic, we can use the resuh of Eq. in Appendix 
C to evaluate 



T T- 

/to 



(21) 



where the factor of 2 accounts for the similar term when the hopping is in the reverse 
direction. Here | ipo) and | Vo) states in the ground manifold with one hole per site, 
and the superscript (2) indicates a result which is second order in t. In using the result in 
Appendix C we must truncate the matrix element so that it remains within this manifold. 
Also, in evaluating this expression it is convenient to use the identity 



'^loo-CjOo-' 



2 + 



cr 



(22) 



whereby we obtain the result 



+ 



b \ 

2Tr{t^{,[l + S(j) 



2 + S(0 • 



Uo + Ejb) 



a 



^00 



+ S(z) - a 



(23) 



denotes the sum of all 



where the traces are over the 2x2 matrices in cr-space and {i ^ j] 
previous terms with i and j interchanged. 

The first term in (^3|), which only involves hopping of a single hole (from site i to j 
and back), is easily shown to be independent of the spins at i and j. [This follows directly 
from the identities of Eq. (pIS]), or more simply from time reversal invariance.] Therefore, 
this term contributes a spin-independent constant, and does not affect the splitting of the 
ground state. Similarly, the terms coming from the factors of 1/2 inside the square brackets 
in the second term also give constants. To order t"^ we have thus arrived at an effective 
magnetic Hamiltonian of the form of Eq. ([12D, with 
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In view of Eq. (T^), this becomes 

n^'\^,j) = ^Tr{(A^i, + tBg . a) (S(j) • a) (A^oo + ^B^oo " ^) (S(0 • ^)}- 



(24) 



(25) 



The symmetry of this form is further discussed in Appendix D, where we show that in fact 
T-C^'^\i,j) is of the isotropic form of Eq. (1). 

B. Contribution of Order t^AW. 



To calculate the contributions of the Coulomb terms of Eq. ([T7|) to the magnetic exchange 
we need to carry out third order perturbation theory. By taking two factors of the hopping 
matrix element we generate terms of order t^. We must include an additional factor of ATic- 
This factor is only relevant in the intermediate state when there are two holes on the same 
site. The relevant matrix element for third-order perturbation theory is written in Eq. ( |(J4D 
of Appendix C. In using this result it is convenient to use the identity of Eq. (p2D. Then we 
obtain the correction to the energy at second order in t including perturbatively the leading 
Coulombic contributions (which we indicate by the superscript "(2,c)"): 



n 



{2,c), 



1 



^CTo-^^g 5 (^0 + Eia) (Uq + Eib 



0-35 



X 



AUss'a2a{i; aObO) + Kss'a2u{i] a060) - MJss'uuii.^ «006) - Kss'au2{i] ^006) 

(26) 



X I i + S(0 ■ a)^^, + S(j) ■ a)^ ^^ + {i^J 



0-10-3 



where we have used the 

property AU„^,„^„'^{i; abb'a') = AU„,„^'^^^{i]baa'b'), K„„,„^^>^{i; abb'a') = k^,^„f^^^{i]baa'b'). 
(In writing the above result we set Ei^Q = Ej^q = for simplicity.) In order to carry out the 
spin summations, we insert here the explicit expressions for AU and K, Eqs. ([18|). This 
leads to 



AUnn' 



Tr{(i + S(z) ■ a) + S(j) • a)x^4 



- Tr{(- + S(z) • a)(mL,o)^mU}Tr{(- + S(j) ■ a)w^^} 
+ [Tr{(^ + S. ■ + S^^') ■ ^)<^'] 



Tr{(- + S(.) ■ a)(mj,,o)tmj,o}Tr{(- + S(j) ■ a^t] 



+ (^ ^ 3] 



(27) 



where 
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xL' = E Tf— T^C(mL) VL,o ^ XZ^J + zY^, ■ a, (28a) 

and 

wii,. = E (ri If uti If ^ ^ol(mL)Vuti^ ^ l^i^J + .Z^^, ■ a, (28b) 

ab y'-^O + -C^mj (tyo + J^ib) 

in which X;^*^, and W^l^, are real scalars and Y^*^, and Z;^*^, are real vectors. It is straight- 
forward to verify that in Eq. (|27| ) the terms which involve two traces, as well as those 
coming from the factors of 1/2 and involving one spin variable, do not contribute to spin- 
dependence in the spin Hamiltonian. This follows by noting that (a) (m^,g)"l'm^,Q and w^'^ 
are proportional to the unit matrix, and (b) one can interchange a and a' in the sums. 

The full effective magnetic Hamiltonian, to order t^, is obtained by combining 'H^'^\i,j), 
Eq. dD, with n^^'^^ij), Eq. (P). The result has the form of Eq. (|), with 

Ji^.j) = 7^Tr{(Poo)(to{,)} - E Af^-'Tr{(xit„0^(xiy^ + 

- E^-'Tr{(xL)'(xit.')' + (xL)'(^-')'>' (29a) 

aa' 



D(^,jO = -7f(Tr{C}TrKV} 



+ ^EA^« 

aa' 

+ 2 ^ 



Tr{x^t,}Tr{(xi^„,) V} - Tr{(x^„,)'}Tr{x^,a}) 
^{xL'}^{(xL)'^} - ^^{{<a')^m<a'^} 



, (29b) 



M(z,jO = 7^(Tr{«®Tr{ti,V} + 



-^EA^. 

aa' 

- 2 E ^a"' 



Tr{xtc?} ® Tr{(xi^,,) V} + Tr{(xi^,,) V} ® Tr{x^^c?} 
Tr{x^,,<?} ® Tr{(x^„,)^a} + Tr{(xi;,) V} ® Tr{xit„a} 



+ (i ^ j] 



+ E Kaa' [Tr{(m*,,o)^ni^o<T} ® Tr{w^t,,<T} + 



(29c) 



One notes that when the contributions of AUaa' and Kaa' are ignored, Eqs. ( ^91) reproduce 
Eq. (i, with D(^,j) = -z[Tr{PoaTr{tj,i,a}/f/o - (z ^ j)], and J = 2Tr{PootoU/t^o. The 
results ( p9D hold for general site symmetry, and to all orders in the spin-orbit coupling. 
They become particularly simple in the special case of tetragonal symmetry, as is discussed 
in Sec. IV. In Eq. ( |29a] ) we see that even when U is a constant and K = 0, J{i,j) does 
depend on A. Moriya's expression for J{i,j) is only correct to zeroth order in A. 
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C. The copper oxygen copper bond 



Here we derive the effective magnetic Hamiltonian of the copper spins for the bond Cu- 
0-Cu. The spin-orbit interaction on the oxygen is much smaller than that on the copper, [ 
PT| and therefore may be neglected. Then the microscopic Hamiltonian (^) is modified as 
follows. First, the kinetic energy now represents hopping between the oxygen and the copper 
ions. That is, in place of Eqs. (0) we now have 

^hop = E^'/* + h- c. , (30) 

iq 

where 

Tq^ =T.T.{tna) y^n^C^a.' , (31) 
an era' 

in which are the destruction (creation) operators for a hole on one of the states 

in) of the gth oxygen, and 

(e)„,=ECK)„,- (32) 

a 

Here describes hopping from the ath orbital on the ith copper ion to the nth orbital 
on the gth oxygen ion. The matrix element (t?^) describes hopping between the copper 
states [see Eq. (|])] which diagonalize 'Hn,. + Ji^^ and the gth oxygen ion. 

Second, we add to the Hamiltonian the on-site single-particle energies and the Coulomb 
interactions on the oxygen. These terms are written in the form 



"^P E ^nP\naPqni7 + „ ^ ^ ^p'^ PqnaPqn'a'Pqn'ij'Pqna 
qna ^ qnn' aa' 



2 ^ ^ ^^nn'PqnaPqn'a'Pqn'cr'Pqna + 2 ^nn'PqnaPqn'cr'Pqna'Pqn'a 
qnn' aa' qnn' uu' 

= n^^^ + An^f^ , (33) 

where Ti^'^ is the first line of this equation and ATi^^ is the second line. Thus, the total 
Hamiltonian for this case is taken to be Tio + V, where 

no = n, + nso + nco + ni^^ (34) 

and V, which we treat perturbatively, is 

V = Htop + AHc + AH^^ . (35) 

In the above, the index q distinguishes between oxygens on the bond along the x- and 
y-directions from the copper ion in question. However, the perturbation expansion gives 
results in the form of contributions summed over all pairs of single bonds between nearest- 
neighboring copper ions i and j. Then, the index q is fixed once the values of i and j are 
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specified, as one sees from Fig. 1. Accordingly, we liencefortfi omit tlie index g, so tliat, for 

We now turn to the perturbation expansion, from which we obtain the magnetic Hamil- 
tonian of the copper spins. It is clear that the lowest order contribution to the effective 
interaction between two copper spins is of order t^. There are two possible channels in this 
order, which we denote by a and h. In channel a, the hole is transferred from one of the 
coppers to the oxygen, then to the second copper, and then back to the first copper via the 
oxygen. Hence in this channel there are two holes on the copper in the intermediate state. 
In channel 6, the hole is transferred from one of the Cu ions to the oxygen, and then a 
second hole is taken from the second copper to the same oxygen. Afterwards the two holes 
return to the coppers, i.e., back to the ground state in which there is one hole on each Cu 
ion. Thus in channel h there are two holes on the oxygen in the intermediate state. When 
the terms coming from the Coulomb interactions At/„„/ and K^n' for the oxygen ions are 
included, then their effect will appear only in channel 6, in which the two holes have a state 
where both are on the oxygen. 

It turns out that for channel a all our previous expressions, derived for the Cu-Cu bond, 
hold with the replacement 

t^. = E ^ • (36) 

We show this explicitly in Appendix E for the process. Similar arguments hold for the 
processes of order t'^AU and t^K, where AU and K are the Coulomb interactions on the 
copper ion. (Note that for this channel AU and K represent the Coulomb interactions on 
the copper). 

It thus remains to investigate the perturbation expansion in channel b. Applying once 
to ipQ the term in the Hamiltonian of Eq. (|30[) , which describes hopping from the copper 
ions to the intervening oxygen ion, one obtains 

I V'l > ={Ti + Tj) clo^c]Q^,CjOa,Cio„ I V'o) 



n aa\a2 



which represents virtual states with energy e^. In the next order, the second hole is put on 
the same oxygen. This leads to 



\/to /To lio rCo 



= EE E (f + 7;) , +,\ + u iQ.,M'o).,.Pn^Ms^^o.,c,o. I V^o) . (38) 

nn' o-ai 0-20-3 t„' t„ T t„' T Up 2 3 1 

In order to return to the ground state two more powers of the hopping are needed. This 
gives 
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„„i ^n' ' ^n' 



Tr{(l + a-S(0)t^Xo} 



Tr{ (l + a- S{j))ti^,ti,,} - Tr{ + a • S(^))tj,„,P„,o(^ + a ■ S(j))Pont^o} 



(39) 



where we used the identity Eq. (0). (We labeled this contribution with a superscript 2 
because even though it is fourth order in the t's, it is really a second-order process in terms 
of a renormalized Cu-Cu hopping interaction. The superscript "b" indicates a contribution 
from channel b.) In a similar way to the arguments given after Eq. (P^D, one can convince 
oneself that the first term in Eq. (|39D as well as the terms coming from the factors of 
1/2 in the second term do not contribute to the spin Hamiltonian. Thus, to order t"^, the 
contribution of channel b is 



^(26) ( 



nn' ^n' ~r tj^' ~r Up 



Next we calculate the effect of the Coulomb terms Af/„„/ and Knn' of Eq. 
end we apply them to the state \ip2b) of Eq. (pSf). The result is 



(40) 
To this 



— (AU + Kj|V^2.) 



EEE(- + - 



nn' ccri o"2 0"3 " 
t t 



cr2Cr \ / (T3Cri 



(41) 



Finally we apply two factors of the hopping which bring the holes back to the ground state. 
This leads to 



(26,c) , 



1 



Cn' ^ + Cn' + Up 



A[/„„,Tr{a ■ S{{)i^,^rti,,a ■ 8(^)^0X0} 



+Kr.n'Tr{a ■ S(^)t^„,t^o}Tr{a^ • SUftiA,} - J^Trja • S{z)tlX'o^ ■ S{j)ti^,t^] . (42) 

Combining Eqs. (^) and (^) we obtain the magnetic interaction arising from channel 
b in the form of Eq. (El) , with 



+ e„/ + Up 

1 \2 fr„„/ 

e 



(1 



e„ + e„' + ?7„)2 



Tl'{^On^n'o(^On'^no) } 



)'^^{^On'^n'o(^Ontno) } 



(43a) 



1x2 1 

+ e„' + Up 



(1- 



Af/„„,' 



-) 



( Tr{Pont^o}Tr{t^„,F„,o^} " Tr{t^„,F„,o}Tr{FoXo'?} 

f;\rr (Tr{Po.-t-o}Tr{t-„P„,oa} - Tr{t^„P„,o}Tr{Fo„,tt,oa} 



, (43b) 
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^ „„i Cfi' 



Cn + Cn' + Up 



Tr{t^„,F„,oa} ® Tr{Pontj.o^} + TrfPont^o^} ® Tr{t^„,tl,oa} 



2Tr{t^„,t^,oa}®Tr{FoJ^„,o'?} 



Tr{t^„P„,o'?} ® Tr{Fo„,t^o<?} - Tr{Po„,t^,o<T} ® Tr{t^„P„,oa}) 



(43c) 



The full magnetic Hamiltonian for the copper spins of the Cu-O-Cu bond is obtained 
by combining the results of Eq. (^3]) for channel b, with those for channel a given by Eqs. 
(p9|), in conjunction with the identification of Eq. (|36D. These results generalize those of 
Refs. O, n!B| and 1221, which were obtained in the absence of the Coulomb terms AU and K. 



IV. TETRAGONAL SYMMETRY 

This section consists of three subsections. In Sec. A, we apply a canonical transforma- 
tion to show that without Coulomb exchange interactions the effective spin Hamiltonian is 
isotropic at all orders of t and A. In Sec. B, we use this theorem to isolate the most impor- 
tant contribution to the anisotropy, namely that involving the Coulomb exchange energy. In 
Sec. C we corroborate our analytical results of perturbation theory by numerical solutions 
for single-bond clusters: Cu-Cu and Cu-O-Cu. 



A. Canonical Transformation 

We start by proving our strongest result, namely that the spin Hamiltonian Tis arising 
from the generic model is isotropic for a wide class of models in the absence of Coulomb 
exchange interactions. In particular, this result holds for a commonly used model of the 
cuprates, in which interionic Coulomb interactions, V and M in Eq. (^, and the Coulomb 
exchange terms K are neglected, hopping is between nearest-neighboring Cu ions, and the 
site symmetry is tetragonal. Strictly speaking, the only use we make of site symmetry is 
that it has to be high enough so that the 3d spatial orbitals which diagonalize the crystal 
field Hamiltonian are ipo{r) ~ — y"^, ipi{r) ~ 3z^ — r^, ipx{^) ~ yz, "ipyi^) ~ xz, and 
ipziT^) ~ xy. Here the z axis coincides with the tetragonal c-axis and the x- and y-axes 
coincide with the nearest-neighbor directions in the plane perpendicular to the c-axis, as 
shown in Fig. 1. These symmetry labels are chosen so that '?/'a(r) transforms (under the 
operations of tetragonal symmetry) like La for a = x,y,z and ipoi'^) and ipi{r) transform 
like scalars. 

An important observation is that the orbital angular momentum operator L has matrix 
elements only between states of specific symmetry. For instance, connects '?/'a:(r) only to 
the states ipoi'^) and ipi{r) and it connects ipy^r) to '^/'^(r) and vice versa. Similar statements 
can be made about the other components of L. We now introduce a transformation in 
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spin space (to pseudospin) such that the spin-orbit interaction is diagonal with respect to 
pseudospin. For that purpose we introduce pseudospin jl as follows 

|a, /i >= i)a,{r) YM^]l,v^v = /i/.|vac) , (44) 

V 

where 0^ is a spin function for spin "up" if 1] = 1/2 and for spin "down" if rj = —1/2, 
and |vac) denotes the vacuum state. Here for a = x, y, z are the Pauli matrices, and 
(To = 0"! = I is the unit matrix. As discussed in Appendix F, the above transformation is 
such that the spin-orbit interaction is diagonal in pseudospin: 

^so = A 5: (^)fLM, , (45) 

k=i,j a,f3,fi 

where W{k,a, (3) is a spin-orbit matrix element. [The transformation of Eq. (^^ should 
not be confused with Eq. (§). The latter involves an exact diagonalization and requires a 
knowledge of all the parameters. In contrast, the transformation of Eq. ( |^ ) is independent 
of the interaction parameters. It merely brings the Hamiltonian into block diagonal form 
in which there are two identical blocks, one for = 1 and one for fi = —1.] Because the 
unitary transformation of Eq. ( |1^ ) does not mix spatial states, it does not affect the form 
of the Coulomb interactions scaled by U. Furthermore, in view of the lattice symmetry 
hopping can only involve holes moving from one site to a neighboring site without changing 
their symmetry. Thus holes in a state ipa on one ion, where a = x,y, or z, can only hop to 
states of the same a on a nearest-neighboring ion. Likewise holes in states V'o(r) or ipii^) 
on one ion can only hop to states V'o(r) or ipi{r) on an adjacent ion. Since states a and (3 
which are connected by hopping must be states of the same symmetry, we have 

'^ij ~ 5Z ^ia,jl3dlafj.djl3^i = ^ ^«a,i/3 [^alr^i/iLr [^/3]p/i/i/3p — 5^ tiaJ/sfiapfjPp ■ (46) 
Q,f3,fi a,l3,fj.,p,T <^,P^P 

In other words, the total Hamiltonian (for K = 0) can be written in the form 

^ = ^kcxfkapfka^i + X! W{k,a, P)fla^fkp^i + UaJ P fLpfj Pt^ 

k Q^iA* k Q,f3,fi i,j,a,l3,fi 

Uka,ka' fkap,fka' p' fka' fi' fkap- (47) 

k a,a' ,p,p' 

Thus this Hamiltonian can be written [ |2^ in terms of the quantities 

Q«/3(^,j)=E/l/./.73M (48) 

which themselves are invariant under rotations in pseudospin space. Therefore Ti. is invariant 
under rotations in pseudospin space. To construct the effective spin Hamiltonian H{i,j) 
involves using degenerate perturbation theory to eliminate the excited states. Accordingly, it 
is clear that the resulting spin Hamiltonian will be rotationally invariant in pseudospin space. 
Since we have defined pseudospin so that in the ground state (in which all the holes are in the 
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state ^o) pseudospin and real spin are identical, it follows that for the tetragonal case with 
no Coulomb exchange interactions the spin Hamiltonian is also rotationally invariant. This 
theorem indicates that even though Tis may include further neighbor two-spin interactions, 
four-spin interactions, etc., it is nevertheless rotationally invariant, so that the spin-wave 
spectrum can not have a gap at zero wave vector in the absence of Coulomb exchange terms. 

For the case of the nearest-neighbor exchange interaction we can make some further 
explicit statements. For tetragonal symmetry, J^u{,h3) must also be diagonal (with its 
principal axes along the tetragonal axes). Thus, T-C{i,j) is an isotropic Heisenberg model. 

The above theorem can be generalized to include the intervening oxygen ions. Here we 
consider a Hubbard model which includes the three 2p spatial orbitals. Now we introduce 
different unitary transformations for oxygen ions on y-directed and x-directed bonds (see 
Fig. 1). For those on ?/-directed bonds we set i^o^r) = \2py), ipxi'''') = |2pz), and '^^(r) = 
|2px)- We then introduce states by 

l«, /^) = V^a(r) J2M^,,r,M^) ' (49) 

similar to Eq. (H) which was used for the Cu d states. We need to examine how the hopping 
and spin-orbit interactions are affected by this transformation. Note that hopping along the 
y-direction can only take place between Cu states like — y"^ and oxygen 2py states. These 
are both associated with symmetry or 1. Likewise an oxygen 2pz orbital can only hop to a 
copper yz state, both of which have symmetry label x. Also an oxygen 2px orbital can only 
hop to a copper xy state, both of which have symmetry label z. Thus with this labeling of 
states, hopping occurs only between states of the same symmetry label and the canonical 
transformation has no effect on the hopping, just as in Eq. (|^). One can verify that the 
spin-orbit interaction on the oxygen ions does conserve pseudospin. Oxygen ions on the 
x-directed bonds are treated analogously. For them we write '?/'o(i") = \2px), V'j/(r) = I'^Pz), 
and ^'^(r) = \2py) and we again use Eq. (^9]). Then, we conclude that the Hamiltonian of 
the entire lattice can be expressed in terms of the quantities Qapihj)- Thus the theorem 
holds with intervening oxygen ions: in the absence of Coulomb exchange, this model gives 
no anisotropy in Tis for a tetragonal lattice. 

The fact that this conclusion is demonstrated to all orders in perturbation theory repre- 
sents an important new result. The low order perturbation result of BS is in accord with this 
theorem. As mentioned there, this conclusion modifies the conventional wisdom that the 
anisotropy in the exchange interaction is trivially related to the anisotropy of the g tensor. 
Finally, we emphasize that this theorem depends crucially on the fact that the eigenstates 
of the crystal field are those of tetragonal site symmetry and that hopping is only between 
nearest neighbors. In addition, the theorem is only valid for Coulomb terms which have the 
Hartree form, i.e., those which only involve two orbitals, as in Eq. (^). More complicated 
Coulomb terms (involving more than two orbitals) and Coulomb exchange terms [involving 
K in Eq. @)\ will generate anisotropy. 
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B. Anisotropy 



Anisotropy can occur via various mechanisms. One such mechanism is to introduce 
Coulomb exchange interactions, as done imphcitly by BS. In the context of the above dis- 
cussion we note that exchange interactions compete with spin-orbit interactions in the fol- 
lowing sense. With only the former interactions the eigenstates of a single ion are states of 
total real spin 1 or 0. With no exchange but with spin-orbit interactions, the eigenstates of 
a single ion are states of total pseudospin 1 or 0. In both cases, our numerical evaluation 
of the energy levels gave singlets and triplets as this argument requires. However, the wave 
functions are different, of course. When both interactions are present, the degeneracies are 
removed because wave functions can not be simultaneous eigenfunctions of both real spin 
and pseudospin. /^From our calculation, treating hopping, spin-orbit. At/, and exchange 
interactions as perturbations, we reach the following conclusions. For the "generic" model, 
anisotropic exchange appears at order t^A^K. An efficient way to perform this calculation 
is to use the hopping matrix elements, t, of Eq. ( |14b|) and work to order t^K, as we did in 
Sec. III. The perturbation expansion yields the general expression, given in Eq. (p7D. In 
Appendix G we analyze this expression for tetragonal symmetry and find that it agrees with 
the result given previously [ [T^ , namely 




+ ei + f/i,^) 



1 1 
— + FT- 



ei + f/o,i J' 

where Lj^^ denotes the orbital angular momentum matrix element, (a|L^|/5), /i assumes the 
values X, y, and z, and the superscript "anis" indicates that we have arbitrarily omitted 
isotropic (i. e. /^-independent) contributions. The same expression is also derived directly 
from perturbation theory in t. A, and K in Appendix H. 

Now we briefiy discuss the implications of the above result. First of all, note that within 
tetragonal symmetry the result does display the expected full anisotropy for a single bond, 
under which J^x = J\\, Jyy = J±, and Jzz are all different. To get biaxiality (Jy 7^ J±) 
requires either toi 7^ or t^x 7^ tyy, viz. Fig. 2. Of course, in tetragonal symmetry single- 
site quantities can not differentiate between the x (||) and y (±) directions. To understand 
why toi 7^ introduces biaxiality, note that toi changes sign when the local x coordinates 
are rotated into the y coordinates. Also note that even in the limit when U is considered 
to be very large, the result still does depend on the hopping between excited levels through 
t^^. Finally, we remark that these expressions differ in several respects from those of BS. 
This point is discussed in Appendix I, where we give the results more explicitly. 

C. With Oxygens 

Turning now to the Cu-O-Cu bond in the tetragonal symmetry, we again discuss sepa- 
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rately the anisotropy resulting from channel a (the two holes occupy the same copper in the 
intermediate state) and that coming from channel h (the two holes are on the oxygen in the 
intermediate state). 

For channel a, we use the transformation of Eq. (PB]). Using Eqs. ( |G5| ) and Eq. ( |52D we 
obtain 

\ / a^ (JO a _ ^ \ / cri cro 



na/3 



where scalars [see Eq. (p5|)]. For tetragonal symmetry, a and /? belong to the same 

symmetry class (e.g., a = P for a = x,y, or z, or a and /3 are or 1). Hence tj'^ becomes 
the unit matrix times a scalar given by 

a/3 n 

in which a and /? refer to tetragonal d states. The calculation can now proceed exactly as 
for the "generic" model described in Appendix G (or H), with the effective hopping matrix 
elements t^/s given by 

2 / / 

^00 = tQp^/ ep^, toi = topjp^i/ ep^, 

txx = 0, tyy = iyp^/^p^J ^ ZZ = ^ zpy l ^Py ^ (53) 

where px-, Py, and pz represent the p states on the oxygen. Thus the contribution of channel 
a to the symmetric anisotropy of the spin Hamiltonian for a bond along the x-direction 
is reproduced by Eq. ([5D|), with the replacements (^). Analogous expressions hold for a 
Cu-O-Cu bond along the y-direction. 

Now let us consider the magnetic anisotropy in channel b. Inspection of Eqs. ( [43b|) and 
( [43c| ) shows that we need to examine the 2x2 matrices to„t„/o. (We omit the site indices 
which are irrelevant for the tetragonal symmetry) . Using Eq. (^) and tetragonal symmetry, 
we write 

(54) 

aP ex 

Therefore, they are proportional to the unit matrix. As a result, there is no contribution to 
the magnetic anisotropy in channel b in order t^. The reason is that in this channel (and 
to this order) the excited states on the copper are not visited at all. Therefore, just as is 
the case for the Cu-Cu bond when those states are ignored [cf. Eq. (|2^) ] the magnetic 
Hamiltonian resulting from this channel is isotropic. 



D. Numerical Study 



We have checked our analytical results of perturbation theory against results (shown in 
Fig. 3) obtained from exact diagonalization for the four lowest levels out of the 190 possible 
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two-hole states for a pair of Cu sites. The relations between the exchange constants and 
the four lowest levels are obtained in Appendix A. 

For our numerical results shown in Fig. 3, we used the values of the parameters listed in 
Table I. [ |2^, |2^, ^ ^ The hopping matrix elements are related to (pda) as 

follows: to,p, = -V^h,p, = ^ipda) and ty^^^ = t^^p^ = {pdn), with (pdn) ^ -i(pdcr).[ |5§ 
The expressions for Ua^j3 and K^^/s in terms of the Racah parameters were taken from Ref. 



29| and are listed in Appendix B. We also checked that the J^J^^*^, shown in Fig. 3, agree to 



within about 10% with those obtained from the full 325 site Hamiltonian for the Cu-O-Cu 
cluster. Very crudely, as A increases above 7 eV, the are inversely proportional to A~'^ 
and are proportional to a linear combination of B and C. (When B = C = 0, our theorem 
indicates that there is no anisotropy in J^^.) Thus our results are not highly sensitive to 
increasing the value of A. As A is decreased below about 6 eV, perturbation theory rapidly 
becomes increasingly inaccurate. Now we discuss briefly the numerical values of the Racah 
parameters. We took the values of B and C from Ref. 31. Then, fixing the value of A is 
equivalent to fixing the value of Uq = A + AB + 3C. Recently proposed values for Uq are 8.8, 



8.8, 9.4, and 10.5 eV from Refs. ^ and |2^, respectively. As a compromise, we took 

Uq = 9.34 or y4 = 7 eV. Our parameters yield an anisotropy in J of order 0.03 meV and, as 
we shall see in the next section, give an out-of-plane gap in the spin-wave spectrum within 
10% of the experimental [ p| value 5 meV. 



V. SPIN WAVE SPECTRUM OF THE EFFECTIVE SPIN HAMILTONIAN IN 

TETRAGONAL SYMMETRY 



Given Eqs. (|I3|) and ( pOD for single bonds, the classical ground state of the effective 
spin Hamiltonian is rotationally invariant in the basal plane. The out-of-plane anisotropy 
axY ~ ^J/Jq is positive (see Fig. 3), and therefore the spins order in that plane, as is well 
established. In the absence of spin wave fluctuations, the in-plane gap is zero. However, 
the classical rotational invariance within the basal plane is broken by the dependence of 
the spin-wave energies on the angle 6 between the staggered magnetization and the crystal 
X-axis. The purpose of this section is to study this anisotropy and show that it leads to a 
nonvanishing in-plane gap in the spin-wave spectrum. [ |l£ 



In order to show this, we start with the following general Hamiltonian for the Cu02 
plane in a tetragonal system 

H,s = J2 H^J , (55) 

<ij> 

where for < ij > along the x-direction, Hij is 

= J|| S^S] + J±SfS] + JzS^S] (56) 
and for < ij > along the y-direction Hij is 

Hi, = J±S^S^ + J|| SfS] + JzS^S] . (57) 
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We will now calculate the spin-wave spectrum of this Hamiltonian and then the first 
quantum correction to the classical ground state energy. We consider the case where the 
spins lie in the xy-plaiae and ordered antiferromagnetically {J\\, J± > Jz > 0). Assuming 
the staggered magnetization moment makes an angle 9 with the positive x-axis, we use the 
following transformation so that spins are parallel to the new 2;-axis: 



/ — sin 6* cos 9 \ 



V 



cos 9 sin 9 

-10 y 



(58) 



Defining sublattice A to have up spins (in the rotated frame) and sublattice B to have down 
spins (in the rotated frame), we have the following bosonic spin representation: 



ai + a'J 



IS 



aj — a- 



(59) 



for sublattice A, and 



s'; = -s + bp, 



(60) 



for sublattice B. For later convenience we consider the case of general spin, although in the 
end we set S = 1/2. Using Eqs. (|S^-^D|) we may write the effective spin Hamiltonian H^s 
given in Eq. (|55|) in momentum space as 



Hes = Eo + AJ^^S aJjOq + bpci + A^ap^^ + B^ap\ + h. c. 



(61) 



where q is summed over the first Brillouin zone of the magnetic reciprocal lattice and 



1 



>/av = -{J\\ + Jl 



Br 



4J. 



av 
1 



4X 



[Jl cos^q^a) + J2 cos(gj^a)] , 
-{.hcos{qxa) + JiCOs{qya)\ 



(62) 



Here is the total number of spins and 



Jl = J|| sin^ 9 + cos^ 9 + Jz , 
J2 = J|| cos^ 9 + J± sin^ 9 + Jz , 
J3 = J|| sin^ 9 + J± cos^ 9 — Jz , 
J A = Jl cos^ 9 + J \ siv? 9 — Jy. 



(63) 



Henceforth we will set the lattice constant a to unity. Note that our conventions imply that 
X]q 1 = A^/2. As one expects the classical ground state energy Eq does not depend on 9 
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and thus we have complete degeneracy with respect to 6. However diagonahzation of the 
Hamiltonian in Eq. ( |6lD leads to the result 

i/eff = ^0 + H I ^+(q)aqaq + ^-(q)^q^q f > (64) 
q J 

where the new ground state energy Eq is now 

K = -2(1 + ^)NJ,^S' + 1 ^ |^+(q) + a;_(q)| (65) 

and thus does depend on 9. This dependence on 9 arises because the zero-point motion 
contribution (which is the sum of spin-wave energies uj^{q) +ci;_(q) over the Brillouin zone) 
depends on 9. The spin- wave energies are 



^+(q) = 4Ja.5^(l- 5^)2-^2 

(q) = 4J,.5y'(l + 5q)2-A2 . (66) 

Note that when Jy = J± = Jz, -Bq is zero and thus we have two degenerate spin modes 
as usual. When Jy, J± and Jz are different, the two modes are no longer degenerate. This 
remains true when q — 0: 



^+(0) = AS^2JM, - Jz) ; 00.(0) = . (67) 

This result shows that we have only one gap in the noninteracting spin wave picture even 
though the ground-state energy is anisotropic and therefore selects [ Q a value of 9. 

In Fig. 4 we plot the noninteracting spin-wave spectrum according to Eq. (|66D along 
different directions in the Brillouin zone. For illustrative purposes we arbitrarily chose values 
of the J's which correspond to much larger anisotropy than we have for the cuprates. An 
indication [ ^ in the spin-wave spectrum that 6 J = J\\ — J± is nonzero is the removal 



of degeneracy between co'+(q) and a;_(q) on the boundary of the Brillouin zone (where 
Qx + Qy = tt). This effect is illustrated in Fig. 4. Even though noninteracting spin-wave 
theory does not lead to two gaps at zero wave vector when Jy ^ J±, one can obtain the 
second gap by calculating the spin-wave spectrum including higher orders in 1/5*. However, 
below we will estimate this in-plane gap without explicitly invoking spin-wave interactions. 
For this purpose we study the quantum zero-point energy (per spin) in detail. It is given 

by 

^^(^) = ^E(^+(q) + ^-(q)) • (68) 



/^From Eqs. (|62D, (|63D, and ( |66D one can write 

5J' 



^^±(q) = 4Jav5' 



if±g) + ih±k)cosi29)(^—^ 



1/2 

(69) 
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where 



4 7 

(yJa.v ~ Jz)C-\ 



2X 



8Jav 
(J_ 

k = — — , 6 J = J|| — J_L 



with C+, C„ given by 



C± = cos{qxa) ± cos(gj,a) . 



(70) 



(71) 



To obtain the leading ^-dependence of the mode energies, we expand uj± (q) up to second 
order in powers of (5J/Jav): 



.,(q)^4J..S(/±,)V^{l + i(^)cos(2»)(^ 



6J\ l/h±k 



J,J 8\f±g 



cos 



'Mm 



(72) 



By using this in Eq. (^) we can obtain the leading 6'-dependence of the quantum zero-point 
energy, 

(sjr 



where the numerical constants are 



SJ 



Co + Ci cos(2^)^ - C2 cos2(2^)- 



(73) 



Co 
C2 



1 f + fc) (/i- A:) 



+ = 



(/i + fc)^ 



(/ + 9)'/' ^ if- gf/^ / ■ 



9- 

(h-k)^ 



\ 



(74) 



Note that coefficients of odd powers of 6 J vanish due to the fact that these terms include 
odd power of C_ which changes sign under Qy while the other expressions are invariant 
under this operation. 

In Fig. 5 we show Ez{0) from Eq. ([73|) and from the exact sum given in Eq. (^) for 
J|l = 1, = 0.9 and J, = 0.8 for which Cq = 0.44 + 0.39 = 0.83 and C2 = 2.95 x 10-^ + 0.7 x 
10~^ ^ 1 X 10~^, where the first and second numbers are the contribution from out-of-plane 
and in-plane modes, respectively. Note that the in-plane mode contributes almost twice as 
much as the out-of-plane mode. The agreement between the exact and approximate results 
is excellent even though we have taken SJ/Jav ~ 0.1. Since in many real systems this ratio 
is extremely small, Eq. (^) should give nearly the exact value. For Jy = J± = = J we 
have 



Co = 0.842 , = 1 X 10" 



(75) 
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Note that the zero-point fluctuation energy favors the staggered magnetization to point 
along a [1,0] direction within the easy [ ^ plane. Experiments [ |38| indicate that this 
may be the case, for YBa2Cu306, where the dipolar energy does not select a value of 6*, [|T^ 
although it is not easy to distinguish the direction of the staggered magnetization in such 
systems. [ For other tetragonal cuprates, the magnetic structure in the ground state 
is determined by the competition between Ez{9) and other anisotropies which result from 
inter-plane interactions. [ |19| 

We are now ready to estimate the in-plane gap due to the anisotropy of quantum zero 
point energy shown in Fig. 5. To do this we assume that the quantum zero-point energy is 
equivalent to an effective Hamiltonian for general S of the form 



H^zv^ = Y.l:nS[sfSf/S'^) . 



(76) 



Since we are interested in the region where q ~ 0, this effective interaction is probably 
adequate to approximate the dependence of Ez{9) on di even when 9 has a slow nonzero 
spatial variation. Note that -ffqzpE is of order S because Ez{0) in Eq. ( pB| ) is of that order. 
By comparing the angular dependence of Eq. ( [73| ) and Eq. ([7B|), one obtains 



7in = 8C2-^ = 44v5i, 



(77) 



Transforming Sx and Sy into the local quantization axis by using Eq. (|58D (with = 0) and 
Eqs. (^9H60|), and keeping only the terms at order of l/S**^ = 1), we find that 

1 



-f^QZPE = 4A^ Javf^in + 4 Jav^in ^1 1 '^I'^i + ^i^' 



b; + h.c. 



(78) 



In momentum space, -ffqzPE is 

i^QZPE = 4A^ Jav^in + 4 Jav^in \ 



aqfl-q + 6q6_q + h.C. 



(79) 



Hence the total Hamiltonian Hiot = Hqzpe + H^s, where H^s is given in Eq. (pT]), is 



tot 



-E'tot + 4Jav5' 



^(l + 5-X0(<aq + &i&q 



n ^ ^ n ^ 



(80) 



The spin-wave energies uj± are given as 



^+(q) 



4Ja,v5 



1 - Aq - fiq 



1 + 2S-'6ir, + Aq 



q / ' 



u;!(q) 



4Jav5 



1 + A„ + B, 



1 + 2S-'6in - Aq + 5q 



(81) 
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To get the in-plane and out-of-plane gaps, we set q = in which case 



u+{q = 0) = AJ,,SJ2[1 - (J,/Jav)][l + 5-15, 



uj^q = 0) = 4J^.Sp[l + (J./Jav)]5-l(5in, (82) 

where = 0) is the in-plane-gap due to the quantum zero point energy that we are 

looking for. Note that to oc ^/S, as was originally found [ ^ in a similar situation where 
the gap is due to quantum zero-point effects. Thus we see that the noninteracting result, 
plotted in Fig. 4, which gives one gapless mode, needs to be modified as we have just done. 

We now give a numerical evaluation of the gaps uj±{q = 0). For that purpose we approx- 
imate the result of Eq. (p2]) as 



u^q = 0) = 45 J2Jav(Jav - Jz) 



uj^q = 0) = 8J^,SJS,JS = 86JJ{2S)C2 ^ 0.8 | J|| - J± | , (83) 



where we used Eq. (|77p. To evaluate uj+{q = 0) we use the experimental value [ ^, ^ 
Jav = 130 meV and take Jav — Jz = {J± — Jz) — iJ± — <^||)2 = 30 /leV from Fig. 3, in which 
case 

o;+(g = 0) = 25(5.9meV) . (84) 

Measurements [ |TT|] show that zero-point fluctuations reduce 25" to about 0.8. Using this 
value, we get uj+{q = 0) = 4.7 meV, which compares favorably with the experimental [ |^ 
value of 5 meV. From the data shown in Fig. 3 we see that the in-plane gap, uj-{q = 0) 
should be about 25 fieV. It would be interesting to observe this via an infra-red absorption 
experiment. Because the theoretical estimate of the frequency range is uncertain, it might 
be useful to locate the mode at high magnetic field and follow it back to zero applied field. 



VI. THE LOWEST SYMMETRY MODEL 



In Sec. IV we showed that for the model of Eq. in tetragonal symmetry the anisotropy 
vanishes in the absence of Coulomb exchange. This was due to the fact that we have only 
hopping between orbitals of the same symmetry. However the theorem breaks down when 
we have nonzero hopping between orbitals of different symmetry. Thus in this section 
we consider a system with lower symmetry to show that we can have anisotropy without 
Coulomb exchange. 

Here we again consider the effective spin Hamiltonian for two copper ions, but now we 
do not assume any particular symmetry. Thus the orbitals localized on the two Cu ions 
which diagonalize Ti^ are no longer the same and will be some arbitrary linear combinations 
of — y'^, 3z^ — r^, xy, yz, zx, respectively. We write these orbitals as 

4<. = E^«i(04' (85) 
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where Ra^pii) is the matrix element of the orthogonal matrix which gives the new states in 
terms of the undistorted d-orbitals for the ith copper ion. 
Within these orbitals the hopping matrix elements now are 

'^iajP = X! Ra,'r{'>')Rl3,r]{j)t^,rj, (86) 

where t^^^ is the usual overlap integral between the undistorted d-orbitals, listed above. 
Similarly, the matrix elements of angular momentum in this new basis are 

(87) 



A. Numerical Study 



We now present our numerical results for the effective spin Hamiltonian when we use 
these new hopping and angular momentum matrices for two arbitrarily chosen matrices 
R(i) and R(j). In Fig. 6 we show the anisotropy (energy differences between triplet states) 
for two different situations; (1) the on-site Coulomb repulsive interactions depend on the 
orbitals {Ua,f3 7^ U) (2) Orbital-independent (constant) Coulomb interactions: Ua,i3 = U . 
As in the tetragonal case (but now with no Coulomb exchange interaction), we have full 
anisotropy for both cases. By fitting the numerical results shown in Fig. 6 we showed that 
anisotropy is proportional to t^\^ for nonconstant U and to t^\^ for constant U. In the next 
section we give an analytic proof that for constant U the anisotropy vanishes up to order 



B. Order Results for Constant U and K = 



With only nn hopping on the square lattice, there are no contributions at order t^. 
The calculations at order will generate two types of contributions: [ ^] one a four-spin 
interaction, the other two-spin interactions either between nearest neighbors or between 
next -nearest neighbors. The first type is generated when a hole hops around a closed loop, 
i. e., from site 1 to 2, the hole which had been earlier on site 2 hops to 3 and so on, until 
the hole from site 4 hops to 1. The second type of interaction is generated both by closed 
loop processes and by various arrangements of four hops involving two or three sites. 

In this paper we are mainly concerned with the evaluation of the nn pair exchange 
interactions. In particular we have concentrated mostly on the anisotropy of these interactins 
due to spin-orbit interactions. To study the contributions to this anisotropy from repeated 
hopping within a single bond to order t'^, we use Eq. (C^) of Appendix C and apply to it 
two more hopping terms, TijTi.Q^TjiTCQ^ , ending at a ground state. After some algebra we 
obtain 
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j) = Tr { (a ■ S,) {^2 i^J^ (a ■ S,) (4) } 



2 

X — 



aft 

1 1 

+ 



(88) 



Combining Eqs. (^) and ([88|), we end up with Eq. (^^), in which both {Ai + iBi ■ (?) 



and {A2 + 262 ■ a) are of the form tpg + Ol'^^)- Since at order we had D2 = 0, Eq. ( p5| ) 
now yields D2 = 0(t^), and thus the energy sphtting of the triplet due to D2 is of order [cf. 
Eqs. O and (jD^)] BI/A1A2 = 0{i^), irrespective of the details of Eq. (|S|). Thus, the nn 
magnetic exchange interaction becomes anisotropic only at order t^, and this is correct to 
all orders in the spin-orbit coupling A and for all lattice symmetries. This result is indeed 
confirmed by our single-bond numerical diagonalization, as we showed in Fig. 6. 

We end this discussion with two comments. First, note that the separation of (Ai+iBi-a) 
into a sum of terms of orders t and was only possible because the sums over a and b in 
Eq. (|88D all appeared within one matrix [which appears between (a ■ Sj) and (a ■ S^)]. This 
would not have been possible if we had contributions of the kind X]aTr{(a-Sj)Tf ((T-Si)T2"}, 
representing interference between different hopping paths. Such contributions arise at order 
i^, and generate further anisotropy. [Without them, the analysis of Appendix D indicates 
that the energy levels would be two singlets and a doublet]. Second, note that the symmetry 



only fo^ 



contained in Eq. (|25|) would persist to all orders, had we ignored excited states, allowing 

VII. DISCUSSION AND CONCLUSIONS 



A. Discussion 

It is clear that the role of spin-orbit interactions in causing anisotropy in the exchange 
interaction is an interesting and subtle one. In particular, there has been much controversy 
concerning the way the exchange interaction J^^ depends on the crystal symmetry and 
under what conditions one expects to find a gap in the spin-wave spectrum. For a long 
time after Moriya's seminal paper it was thought that one could neglect M in Eq. (^ 
and that the spin Hamiltonian for a single bond would be anisotropic if the Dzyaloshinskii 
vector D were nonvanishing. It was then observed by Kaplan [ |l^ and by SEA, [ |^ that 



although M is of order and D is of order A, one must nevertheless keep both terms when 
discussing the anisotropy or the gap in the spin-wave spectrum. Two other conclusions of 
these authors were 1) the single-bond spin Hamiltonian 'Hsihj) "was rotationally invariant 
and 2) the overall anisotropy of the Cu-0 plane resulted from a frustration between bonds 
with different values of D. In view of the results of the present paper we are in a position 
to state clearly the conditions under which the first conclusion is valid. In particular, our 
results show that rotational invariance (when the Coulomb exchange K is zero) of the single- 
bond spin Hamiltonian to all orders in t is only to be expected when hopping between 
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excited states is ignored, as the SEA argument does implicitly. We remark that since the 
spin-orbit interaction involves coupling to excited orbital states, it is only nonnegligible 
when the energies of the excited states involved are finite. This being the case, strictly 
speaking, it is not totally consistent to neglect hopping between such states, especially since 
the associated hopping matrix elements are comparable to those involving hopping to or 
from the orbital ground state. Nonetheless, as we have seen, the departures from the SEA 
rotational invariance theorem (due to hopping between excited states) are numerically quite 
small in most cases of physical interest. In fact, for the case of constant U considered by 
SEA, the deviations from rotational invariance only enter at order t^. 

In Moriya's original work to order t"^ it was correct to ignore hopping between excited 
states because he considered the case when U was a constant. In this case, as our results 
in Sec. VI show, rotational invariance only breaks down at order t^, because one has to 
go to that high order for hopping between excited states to come into play. When U is 
nonconstant, hopping between excited states leads to anisotropy in 'Hsi.hj) at order t^, as 
our results in Sec. VI demonstrate. These results thus represent a generalization of those 
by Moriya and by most of the literature which followed him and assumed constant U. 

lYiom. this discussion one might now conclude that spin-orbit interactions would lead to 
anisotropy for the Cu-0 planes as long as one includes hopping between excited Cu states. 
However, the theorem of Sec. IV shows that for the special case when the Cu sites have 
tetragonal symmetry, the generic model of Eq. with only t. A, and U nonzero does not 
yield nonzero anisotropy. The same result also applies to the "real" model including oxygen 
ions. This theorem explains why most previous calculations give no anisotropy for tetragonal 
site symmetry and it emphasizes the importance of including Coulomb exchange terms, K. 
It is then clear why the exchange anisotropy is so small, especially (as noted by BS) when 
compared to the anisotropy in the g tensor. We thus find that for each bond the exchange 
interaction has biaxial anisotropy (J||, J_l, and Jz are all different), where the anisotropy in 
J is of order t1.^^_Q^\^ K , or more correctly, it is of order ^cu-Cu-^^-^ '^^ ^Cu-Cu-^^^' where 
B and C are the Racah parameters which represent deviations from the simple constant U 
Hartree term. 

Even though the single-bond exchange has biaxial anisotropy, the classical ground state 
energy, because it is averaged over bonds along [1,0,0] and [0,1,0], does not select an orienta- 
tion of the staggered magnetization within the easy plane. As we have shown, the anisotropy 
within the easy plane results from quantum zero-point fluctuations. In summary, a complete 
discussion of the anisotropy of the Cu-0 planes requires an interesting study of several novel 
symmetries and the way they are broken by fluctuations. 

B. Conclusions 

We may summarize our conclusions as follows. 

(1) For tetragonal site symmetry, with only Hartree-like direct Coulomb terms, the 
effective spin Hamiltonian is isotropic at any order in the parameters t and A. Inclusion of 
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Coulomb exchange breaks this degeneracy at order t^^^Q^X^K for our generic model and at 
order tQ^_QX'^K for the cuprate system with an oxygen ion between the copper ions. 

(2) Since the easy-plane anisotropy (observed via the "out-of-plane" spin-wave gap 
at zero wave vector) has comparable magnitudes for many orthorhombic and tetragonal 
cuprates, it can not depend significantly on the orthorhombic distortion. Our result, Eq. 
(iO)' yields a biaxial anisotropy in the exchange interaction of order tQ^_QX'^K which can 
explain the observed [ ^ out-of-plane spin- wave gap. 

(3) In the tetragonal case, with the exchange interactions having biaxial anisotropy 
given by Eq. (|50|) , the ground state does not depend on the orientation of the staggered 
magnetization within the easy plane. [As shown in Ref. |19|, this remains true when dipolar 
interactions are included.] However, as we show, [ |19| quantum zero-point fluctuations cause 
an anisotropy within the easy plane which leads to ordering of the spins along the (1,0) 
axes , as indeed was claimed to be observed in YBa2Cu306.[ A rough estimate yields a 
resulting "in-plane" spin-wave gap of about 25 /xeV. An experimental measurement of this 
gap would be very desirable. 

(4) In real crystals, the three-dimensional ordering of the spins is determined by a com- 
petition between the anisotropies treated in the present paper and several other mechanisms. 



such as interplane hopping and interactions, as discussed recently in Refs. |T9|, and |44 . 
These may also affect the estimate of the in-plane gap given in conclusion (3). 

(5) For sufficiently low symmetry and without exchange interactions, the rotational in- 
variance of the single-bond Hamiltonian is broken at order t^X^ for constant U. For arbitrary 
Ua,f3 and sufficiently low symmetry, the single-bond Hamiltonian is not rotationally invariant 
even at order t^A^. 

(6) We have given results for arbitrary symmetry for the effective spin Hamiltonian at 
order including, for the first time, the effects of realistic Coulomb interactions. These 
expressions are valid for the orthorhombic phases of La2Cu04. 

(7) In view of the controversies in the literature concerning the results which include 
spin-orbit interactions we have implemented several checks of our perturbative results. First 
of all, we compared the results given in Eqs. ( pQ]) and (^) with expressions obtained by 
treating both the hopping and the spin-orbit interactions as perturbations. In addition, we 
subjected our analytic results for the tetragonal symmetry case to numerical verification as 
follows. We diagonalized exactly the Hamiltonian within the basis of two holes on either 
a Cu-Cu cluster or a Cu-O-Cu cluster. Then we compared the splittings of the ground 
state manifold (in this case, the lowest four states) with those predicted on the basis of 
our analytic evaluation of the perturbative contributions to the spin Hamiltonian. This 
comparison (see Fig. 3) was made with small enough values of the perturbative parameters 
that we can easily check how the results depend on the parameters. 
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APPENDIX A: EXCHANGE AND 4 LOWEST ENERGY LEVELS 



In this Appendix, we explain what we mean by anisotropy and how we study it by iden- 
tifying the eigenvalues of Ti obtained numerically with those of a general spin Hamiltonian 
such as that given in Eq. (|I^). 

The most general effective spin Hamiltonian for a single bond can be written as 



^ Jll Jl2 Jl3 ^ 

J2I J22 J23 
\ J3I J32 J33 J 



(Al) 



where the matrix J is a 3 x 3 real matrix. There exist two transformations Ri and R2 
which transform J into diagonal form if we rotate the spins: Si = RiS'^^ and S2 = R2S2. To 
obtain Ri and R2 we first obtain the orthogonal matrix O which diagonalizes JJ*, where 
the subscript "t" indicates transpose: 



0*fJJ*)0 = 



(A2) 



where is a diagonal matrix with nonegative entries. For simplicity we assume that all its 
entries are actually positive. Then we define J and to be the corresponding diagonal 
matrices with positive entries. 
Then we set 



Ri = O 



Ro = J*OJ 



-1 



(A3) 



where a = DetJ/|DetJ|. In terms of the transformed spins the Hamiltonian is 



f Jx \ 



Jy 







I. 



S' 



2 5 



(A4) 



where all the J^'s are positive. Obviously, a further rotation could be made to change the 
sign of any two components of S'^ (or S2). So the energy level scheme must be invariant 
under such a change of signs. It also has to be invariant under permutations of the J^. We 
find the four energy levels to be 



A. 



E'o - a{J^ + Jy + J,)/4 



a 



X, y, z. 



(A5) 



This set of energies has the proper invariance under change of the signs of any two J^s. 
For the case of arbitrary low symmetry, we did not try to identify the principal axes, but 
tabulated anisotropics, defined to be Ai — A2 and Ai — A3, where Ai < A2 < A3 are the three 
eigenvalues of the set {Aq,}. 

For bond symmetry, J must be diagonal, so that 
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Ti-cR = -E'o + X! JaSl^aS2,a ■ (A6) 
a 

In this case the eigenvalues are given by Eq. ([A5|) with cr = 1 and the Ja have whatever 
signs they have in Eq. ( |A6| ) . The identification of the J's from the set of eigenvalues of Eq. 
( |A5| ) is not unique, because either permuting the J's or changing two of their signs leaves 
the set of eigenvalues invariant. Thus, identification of the J's with coordinate directions 
requires consideration of the eigenfunctions. For this purpose we write them explicitly: 

i^s = \n -iT)/v^, i^z = \n +iT)/v^, 

V'x = I TT - U)/V2 , = I TT + U)/V2 . (A7) 

Then the eigenfunctions are distinguished by their expectation values: 

IpS ■ {Sl,zS2,z) = {Sl,xS2,x) = {Sl,yS2,y) = "1/4 

: -{Sl,zS2,z) = {S,,xS2,x) = {Sl,yS2,y) = 1/4 , (A8) 

and so forth for the other ipa- Having identified which wave functions (coming out of the 
diagonalization of the 190 x 190 matrix) are which, one can easily deduce the values of the 
Ja- For instance 

Jz = -Xs - A, + A^. + Xy . (A9) 



APPENDIX B: COULOMB INTERACTION PARAMETERS IN TERMS OF THE 

RACAH COEFFICIENTS 



Here we list the Coulomb interaction parameters for the tetragonal symmetry crystal field 
states for a d* configuration in terms of the Racah parameters. [ ^ Here AUajj = —2Kai3, 
where Uaa' = Uq + AUaa', with Uq = A + AB + 3C In terms of the triplet and singlet 
energies given in Ref . [ ^ one has = U + K and U* 



U K. 



K 





dx'^ — y'^ 


dsz'^-r'^ 


dxy 




dyz dzx 


dx^ — y^ 





AB + C 


C 




3B + C 3B + C 




AB + C 





AB + 


C 


B+C B+C 


dxy 


C 


AB + C 







3B + C 3B + C 


dyz 


3B + C 


B + C 


W + 


C 


3B + C 


dzx 


3B + C 


B + C 


3B + 


c 


3B + C 



(Bl) 



In the numerical calculations we used (see Table 1) A = 7.00, B = 0.15, and C = 0.58, so 
that Uo = 9.3A, all in eV. 
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APPENDIX C: MATRIX ELEMENTS NEEDED FOR PERTURBATION 

THEORY 



Here we record some of the matrix elements needed to implement perturbation theory. 
If iV'o) represents any state in the ground state manifold (having one hole per site), then we 
may write 



E 



0"10'20"3 



m 

b ^jb + Uo 



g'30'1 J J 



(CI) 



where Tio and Tij are defined in Eqs. ([19D and ( |14b|) , respectively, and where we set the 
ground state energy of Tio to zero. For results to second and fourth order in t with no 
Coulombic perturbations we need to generate the following matrix element: 



'^2) = Tij I ipi) 

= E 



be 



\Ko) (^ci) 4ccr4Cj0cr2 i^bo) (^co) ^^Iccri^^lbas 



(C2) 



We also need 



^3) = :^AHc:^Tj, \ ^0) 
/to /to 



E 



0-30-1 



(j; 010206) + Kss'a2CTz{j\ 010206) 



0-10-20-3 
ss'a-ia2b 



{Ejh + Uo){Eja^ + Eja2 + Uo) \ 

030-2 

(j; O1O26O) - Kss'a3a2{j', 010260) 

^jais^ja2s' ( '^i0o2Ci0o-i I ^0) . (C3) 

Finally, to get the energy at order t^ATic we need 
I ^4) = Tij I ijjs) 



E < 

CTIO20-30-4 

ss'aia2b 



{Ejb + Uo){Eja-^ + i?ja2 + f/o 



At/s^'^2fT3(j; aia206) + Kss'a2<T3{j\ aiO206) - Af/^s'a302(j; 010260) - Kss'a3a2{j; 010260) 



'v \ t t r 

Oai J^^^Cj0o-^C^-0s'fa2,0 



(^0a2)^^^,4oo4c]os'^ai,0 j ^Cj0o2Cj0oi | V^o) • 



(C4) 



In order to make sure this matrix element connects to the ground state, we had to insert 
the factors 6a2,o and Sa^^. 
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APPENDIX D: SYMMETRY OF THE MAGNETIC HAMILTONIAN 



In this Appendix we analyze the eigenvalue spectrum of a system of two spins 1/2 with 
coupling which is arbitrary except that, for simplicity, we consider the isotropic interaction 
to be dominant. In the presence of antisymmetric exchange interactions one can always put 
the Hamiltonian into the following canonical form: 

H(z, j) = a'S(z) • S(j) + ■ (S(z) X S(j)) + yS(0 • M ■ S(j) , (Dl) 

where fi is a unit vector specifying the orientation of the Dzyaloshinskii vector and M is 
a symmetric matrix. Here we will show that the eigenvalue spectrum of this Hamiltonian 
consists of a singlet and a triplet, if and only if the matrix M is such that 7i(z, j) can written 
in the form 

n{z,j) = aS{^ ■ S(j) + (3n ■ (8(2) x S(j)) + l{S{z) ■n){n- S(j)) , (D2) 

where the coefficients obey the relation 

7 = -a + a I 1 + (/5/a)2 |i/2 ^ ^3) 

for finite a/f3. (The reason for phrasing the condition in terms of Eq. ( p2| ) rather than Eq. 
( [Pl| ) is that the former, unlike the latter, is a unique representation.) In this case, as we 
shall see, the spins can be rotated (about the same axis, but through opposite angles) so 
that in terms of the rotated spins the Hamiltonian looks isotropic. This result shows that 
for this relation between the parameters the Hamiltonian is rotationally invariant, even if it 
is not isotropic. (By isotropic, we mean /? = 7 = 0.) 

As seen in the text, many of the perturbative results have the form 

n{z, 3) = -^Trj (Ai + zBi ■ a) (S(j) • a) [A^ + iB^ ■ a) ■<?)}. (D4) 

Defining the vectors D and D2, 

D = AiB2-A2Bi, D2 = ABa + AaBi, (D5) 

this becomes 

- Ima,-,.{ (1 - . ,) (SO) + .E±|i . ,) (s(,) . s) } 
^ '1 + ^k^¥(i) ■ s(') + -tVd • SO) X s(i) 



+ 



t/o (2^1^2)2^ A1A2 

1 



2(^1^2 



S(j) ■ d) {S{^) ■ d) - (S(j) ■ D2) (S(0 ■ D2 



(D6) 



This is clearly of the general form (^, with the Dzyaloshinskii vector 2D/f/o and the sym- 
metric anisotropy matrix 2(D D — D2 ® Y)2) /{UoAiA2). The most general form for M 
would involve introducing a third linearly independent vector D3. 
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We now show that the eigenvalues of the Hamiltonian of Eq. (|D6|) are a singlet and a 
triplet if and only if D2 vanishes. To see this we study its eigenvalue equation, which, after 
some algebra, can be cast into the form 

(A - xf [(A - x)(A + 3x) - 4 I V2 = , (D7) 

where 

4 , . D'^-Dl^ 4 /„ 1 



It is clear that a triplet occurs if and only if V2 = 0, which, in turn, happens if and only 
if D2 = 0. Q. E. D. Furthermore, we see that in the presence of nonzero D2, the triplet is 
split into a doublet and a singlet. To remove all degeneracy it is necessary to introduce a 
third vector D3. 

We make some further remarks about the case when D2 = 0. One can easily verify that 
the conditions of Eq. (P3|) and (^3) are equivalent to requiring that D2 in Eq. (P6|) vanish. 
We further show now that when D2 = 0, the Hamiltonian is rotationally invariant. For this 
purpose note that 

1+*^:^ ^e^'H d = ^, tan^ = ^;^. (D9) 
2A1A2 cose ' \Dy 2A1A2 ^ ' 

The Hamiltonian ( p6|) with D2 = then becomes 

^aMi + — ^)S'(0 ■ S'(j) , (DIO) 



f/o ' {2A^A2r' 
in terms of rotated variables (equivalent to those of SEA [0): 

a . S(j) = e^i'^-V ■ S'(j>-'i'^-^, a ■ S(«) = e-'i''-^S ■ S'(^)e'i'^-^. (Dll) 

When D2 is finite, the triplet splits into a singlet and a doublet. One may ask whether 
it is possible to perform rotations of the spins such that the antisymmetric Dzyaloshinskii 
term will be eliminated and the Hamiltonian will contain only the symmetric anisotropy. 
This is in general not the case. Returning to the Hamiltonian (p^), we put 



a . S(j) = e'^^-'^a ■ S'(j)e-^t^-, a ■ S(i) = e-'^^-'^a ■ S'(i)e*t^-, (D12) 

where the unit vector a and the angle a are yet to be determined. The Hamiltonian then 
takes the form ( P4|) , with the replacements S{i) S'(z), S(j) S'(j), and Ai — > A'^, 
Bi B'i, (i = 1,2) with 

A'l = cosa =F sina(Bi ■ a), B'l = Bi ± Aiasin a — 2 sin^ — (a • Bi)a. (D13) 



36 



The condition that the Dzyaloshinskii term vanish is therefore [cf. Eqs. 



and era 



A 1 cos a — sin a(Bi ■ a) (b 



a(a ■ Bo 



A2 cos a + sin a(B2 • a) yBi — a(a ■ Bi) 

(D14) 



= a sin2a(y4iy42 + (a - Bi)(a • B2)) + cos2a((A2Bi - A1B2) • a 

Since the vector on the left-hand-side is orthogonal to a, Eq. ( pi4D yields 

(A1B2- A2B1) - a 



Ai cos a — (a ■ Bi) sin a 



Bn 



tan 2a 
a(a B2 



AiA2 + (a-Bi)(a-B2)' 
A2 cos a + (a ■ B2) sin a (Bi — a(a • Bi 



(D15) 



It can be shown that these two equations can be satisfied only when the vectors Bi and B2 
are parallel. Hence the antisymmetric anisotropy can be eliminated from the Hamiltonian 
only for specific configurations. Moreover, the criterion for complete rotational invariance 
of H{i,j) is that D2 = A1B2 + A2B1 = 0, or equivalently, that B2/v42 = —Bi/Ai, which is 
equivalent to the condition that 



{A,+^B,.ay = ^^{A2 



iB 



2 ■ o" 



(D16) 



Returning to Eq. (|2^ ) we note that (tgo)^ 



too, therefore that 'H^'^\i,j) is indeed 
rotationally invariant. This represents an alternative proof for the SEA result, [ which 
holds to order t^, to all orders in the spin-orbit coupling A and for all site symmetries 
providing Uajs = U and K = 0. 



APPENDIX E: PERTURBATION THEORY INCLUDING OXYGEN ORBITALS 



Here we show that the perturbation theory results for the Cu-O-Cu bond, through the 
intermediate state in which the two holes are on the copper ion (channel a) are obtained 
from those of the Cu-Cu bond, with the replacement (|36D . 

We start from the state in Eq. (^). As explained in the text, the index q that labels 
oxygen ions on the bonds along x and along y may be omitted for simplicity. Applying again 
the hopping Hamiltonian yields 

/to V /To rCo / 



EE 



02 0"3 



( ^nO ) ( ^an ) 



£3^2 „t „t 



0-30-2 t J 



CjOaiCiOa \ ^0 > 



(El) 
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where we have written the energy denominators exphcitly. We now concentrate on the terms 
that will eventually contribute to the spin Hamiltonian. To order t^, these are obtained by 
applying two more factors of the hopping, that bring the holes back to the ground state. 
The result is 



T 



1 



EE E 

nn' o-aia2 (73(T4cr5 



\ / (T2*^l ^ ^ ^ (T4(7 \ / (75(74 ^ 

from which it is clear that using Eq. (^) one arrives at Eq. (0). 



CjOaiCiOa \ V'O > 



(E2) 



APPENDIX F: DIAGONALIZATION OF THE SPIN ORBIT INTERACTION 

In this Appendix we show that the transformation of Eq. does indeed make the 

spin-orbit interaction diagonal in pseudospin space. We wish to show that 

{(3,u\L-a\a,fi) (Fl) 

vanishes unless /i = z/ in which case it is independent of fi. First of all, note that L has 
zero diagonal matrix elements, i. e. the above matrix element vanishes when a = (3. There 
are now three cases to consider: i) a = and /3 = 1 or a = 1 and /? = 0; ii) a = 0, 1 and 
(3 = x,y, z] and iii) a ^ (3 but both are x, or z. In case i the matrix element of L is again 
zero, so this case is as desired. In case ii with a = (a = 1 is similar) we express the above 
matrix element as 

$:(^^(r)|L,|^o(r))(a;,) = (^^(r)|L^|^o(r))5^,. . (F2) 

In the last step we used the fact that the orbital matrix element is only nonzero when 7 = /3. 
So case ii is as desired. In case iii we write the matrix element as 

= E (V'/3(r)|^7l^a(r))(ac,) Up) (a^) = (?/^^(r)|L^|'i/'„(r))ze/3^a(5/,,^ , (F3) 

where eap-y is the totally antisymmetric tensor. In the last equality we used the fact that a, 
j3 and 7 are Cartesian indices which are all different. Thus all the types of matrix elements 
are diagonal and independent of pseudospin, as asserted. 
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APPENDIX G: RESULTS FOR TETRAGONAL SYMMETRY 



The only nonzero matrix elements of the angular momentum within the manifold of 
normalized tetragonal d states, | >= d^2_y2, \ 1 >= d^z2^j,2, \ z >= d^y, \ x >= dyz, and 
\y>= 4^ are 

(x^ - y'^\L^\yz) = {zx\L^\xy) = -{xy\L^\zx) = -{yz\L^\x^ - y"^) = {xy\Ly\yz) 

= {x^ — y^\Ly\zx) = —{yz\Ly\xy) = —{zx\Ly\x'^ — y"^) = {yz\Lz\zx) = —{zx\Lz\xy) = i , (Gl) 

{3z'^ - r^\L^\yz) = -{yz\L^\3z'^ - r^) = -{3z^ - r^\Ly\zx) = {zx\Ly\2>z^ - r^) = iVs , (G2) 

{xy\L,\x^ - y^) = -{x^ - y^\L,\xy) = 2i . (G3) 
/^From Appendix F and the pseudospin transformation it follows that 

{uj{a,(3))^^^ = u{a,(3){a^ap)^^^^ , (G4) 

where u was defined in Eq. (^), and uj{a, (3) is a scalar. Turning now to the diagonalization 
of the single-particle, single-site Hamiltonian + Ti^o of the Cu-Cu bond, one finds that 
this can be accomplished by putting 

= CT„maa , (G5) 

where scalars which are determined by Eq. (|TT|) , which now becomes 

The solution of this equation, to second order in A, reads 

(G7) 

The hopping matrix elements are diagonal in the tetragonal symmetry, except for the states 
I > and I 1 >, and are independent of the site indices i and j. It follows that t^l^ are 
scalars independent of the site indices as well, and that the hopping is not accompanied by 
a pseudospin flip. This implies that x^a [Eq. (|28a|) ] is diagonal in that space, and that Xqq,/ 
is independent of i and j. Thus the vector D^, Eq. ( p9b| ), vanishes and the matrix M(i, j), 
Eq. ( p9c| ), is given solely by the K^a' terms, and is independent of i and j. 

iFiom. these arguments and using Eqs. (pS]), one obtains that the symmetric matrix 
M(2, j) (which is diagonal in tetragonal symmetry) is given by 
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M(2, j) = 2 5]i^,„,(Tr{ml,om,oC?} ® Tr{w^t."?} 



= 2 II f/T Ur/ ^ (Tr{mJ,,om«o^} O Tr{m],^m«/bCT} 

- ^Tr{mJ^^m„/o'?} ® Tr{m^,om„fe(T} - ^Tr{mJ^,om„b(j} O Tr{m]^^m^/oa}) . (G8) 
Performing the sums over a and a', this yields 

+ Ky^{mlQm*y^ - m*yQrhl^){rhyQrfi^b - rfi^Qrhyb) , (G9) 

with analogous expressions for the yy and zz entries of M. The next step is to write rhaa 
in terms of the spin-orbit matrix elements, Eq. ( p7| ). In doing this we keep in mind that 
both states a and b cannot be the ground state 0, as they refer to intermediate states of the 
perturbation theory. Therefore, it is sufficient to retain for the coefficient of Kq^ the terms 

-S,A, - X6„,^ + A4.5^ - A^ /'^'t"'"-"' . (GIO) 

— — [^x — ^b){^x — ^a) 

Similarly, the coefficient of Kix is 

-A OiaOib . (Gil) 



The leading order of the coefficient of Kyz is of order X^6yb6za, etc. But then toztyo will be 
proportional to too. Therefore, to order A^, the terms arising from Ky^ do not contribute. 
Collecting terms we find the contribution to the anisotropic exchange as 



•^xx V'lj) = ~O-f^0x 77 — ; r^-^Ox^xO + A77 — ; > I toxtaO toatxO 



2Y- 1 1 u;(a,x)^(x,6) ~ ~ 

^ Uo + eaUo + €b{ex - eb){ex - Ga) 



- oKi^A — — ■ -^101*01 n , l^-L^j 



where we have retained terms up to order A^. Finally we write, using Eqs. ( |14c|) and (|G8| ) 

iox = A— ^-^ — ^too ~ t^x) — Xtoi—^ — ^ , (G13) 

~ ^x 

and put a and b in the sums of ( |G12| ) equal to 1, with toi = ^01 • (These are the only possible 
contributions up to order A^). It then follows that 
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= -8A2 



Qx 



1 



a;(0, x) 



(^00 ^ ^xx) 



uj{l, x) 
Uo + ei 



01 



uj{x, 0) 



f/o + ei 



10" 



(G14) 



This result reproduces that of Ref. |1^. This resuh differs shghtly from that given in Ref. 
18] and Eq. (0) in the text, in that the denominators include just the constant part Uq of 
the Coulomb interactions. The corresponding expression in Ref. ^ [and Eq. (^0])] includes 
instead Uq^^, Uo,i, and f/i,^ and thus represents an expansion, described in Appendix H, in 
K but not in AU. 



APPENDIX H: PERTURBATION THEORY FOR THE TETRAGONAL CASE 



In this Appendix we give an alternative derivation of Eq. (^Oj) based on conventional 
perturbation theory in which we treat hopping, Hhop, spin-orbit, TLso, and the Coulomb 
exchange interactions, TYex, as perturbations. For the "generic" model, anisotropic exchange 
appears at order TC^^pH^^He^- To perform this calculation we therefore need to work to fifth 
order perturbation theory and will arbitrarily omit contributions to the isotropic exchange. 
Of course, Tiex can only exist when there are two holes on the same ion, so the five perturba- 
tions must be arranged so that Tihop and Hex occur in the order 'Hhop'^ex'^hop- In principle 
there are ten ways to insert the two factors of Hso- But some study shows that only if the 
two powers of Hso are separated by Hex does the result lead to anisotropy. So the relevant 
fifth order terms in the effective Hamiltonian are 



H,so ^H.hop ~l~ ^hop ^^^^ 



H 



(HI) 



where S is the appropriate energy denominator. If we write the spin-orbit perturbation as 



a ^holos.h 



(H2) 



then it is easy to see that there are no cross terms, i. e. terms involving V^Va with a ^ f5. 
In addition, hopping from site i to site j and back will give the same result as the reverse 
process. So if hopping from site j to site i is denoted T^, then we may write 



Hi^,J)=2Y: 



V —T - + T-—V 



(H3) 



where the operator Qq, that we need to evaluate is simply 



(H4) 
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There are two channels to be considered for the intermediate state in which Hex operates. 
Channel "0" is one in which site j has orbitals |0) and \a) occupied, whereas channel "1" is 
one in which site j has orbitals |1) and \a) occupied. Then we may define 



Q 



(7) 



V —T - + T--—V 



<o,.'4.o,.'|0) ' 



where 7 = or 1. Then 



(H5) 



(H6) 



a 7 



where 



'^ex,a = -o^«7 [2^2: + • ct] 



(H7) 



Here I is the identity operator and cr-a denotes the sum over direct products, J2a (^aO'a- Also, 
each matrix [Q*^'^) or TY^^^] is a matrix in the direct product of the two spin variables. Any 
operator in this space can be written as a linear combination of direct product operators. 
We define AB via 



\AlS\(j-^Y);a' ,ri' — ^(T,iT'-^r/,T;' ■ 

Explicit calculation of the processes shown in Fig. 7 shows that 



(H8) 



(0) 



X{a\La\ 0)t 
2e. 



f/00 
A(a|L„|l)toi 



+ 



00 



aa r r -1 



(H9) 



To write this in operator form, note that | XX + a ■ a 
and is zero otherwise. Thus 

Qi°)] = C, [la^] + [la^] [II + a ■ a] 
= {C^ + C2)[I<J^]-C2[Ia^][0] , 

where [O] = [II - a ■ a]/ 2 and 

^ X f too{a\L^\0) too{a\L^\0) \ 
' n C/00 J 



is unity if cr = 77' and r) — a' 

<T,'n;a' ,r]' 



(HIO) 



(Hll) 



C2 = 



A [ too(Q;|L^|0) _^ taa(Q:|iva|0) 

2 [ UoQ ea 



^ toi{a\La\l) 
(ei + Uio) 



(H12) 



Also 
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(1) 



\toi{a\L^\0) ( 1 



1 



(H13) 



so that 



Q 



Xtoi{a\L^\0) 



(H14) 



Thus we have the result 



X [CiII + C2II - C2O] 



[C{II + CllI - CIO] [Xa^] [Il+a- a] [la^ 



(e„ + ei + f/i«)2 

To simphfy the above result we use the identity for Pauli matrices, 
[XcTq,] [XX + a ■ a] \I(Ja\ = [XX + 2o"q(Tq, — a ■ a] . 



(H15) 



(H16) 



^From the form of Eq. ( [H15|) we see that all the anisotropic contributions come from the 
term aaC^a in Eq- ( |H16| ). Keeping only such terms we have 



{2K 
[C*XX + C*XX - CIO] [a^a^] \C{IZ + C2XX - C^O] 



2K,^\Csmaaaa] 
(ea + ei + f/i„)2 



(H17) 



The terms involving the operator O give only isotropic terms. This can be seen by using 
the equality 



2 [aaCTa] [O] = [(TaCra] [[II - a ■ o] = [u ■ G - II] 



(H18) 



which is isotropic. Thus the anisotropic exchange terms are correctly given by 7{{i,j) 
(1/4) J^T^,^, = J^fS,{z)S,{j), with 



rams 



I ^0,fi I '^O,!-'^!,, 

(e^ + ei + f/i,^) 



1 

— + 



1 2 



+ 



ei + f/( 



0,1 



(H19) 



where L'^g denotes the orbital angular momentum matrix element, 
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APPENDIX I: EXPRESSIONS FOR J, 



ANIS 



In this Appendix we give expressions for J^^^^ [see Eq. (pO])] assuming the relations for 
the hopping matrix elements implied by Eq. (|53| ) and the relations involving [pda) and 
(pdn) listed in Sec. IV. D: 



txX — , tyy — tzz " ^Oo/S, ^01 " ~^Oo/ . 



(11) 



Also, we use the identifications of the Racah coefficients given in Appendix B. Thereby we 
obtain 



B + C 



+ 



3(e,. + ei + A + 25 + C)2 (e, + A - 25 + C)^ 



1 



1 

^e^ ^ ei + A-AB + C 



(12) 



rams 

yy 



-2XH, 



2+2 



B + C 



3(e, + ei + A + 25 + C)2 
35 + C 



1 1 

ei + A-AB + C 



+ 



% + A-2B + Cf 



2 

3e^ 



ei + A - 45 + C 



(13) 



Jzz — — 8A t 



2.2 



45 + C 



00 



+ 



3(e, + ei + A-45 + C)^ 
C 



+ 



1 



ei + A-45 + C 



;e, + A + 45 + cy 



n 2" 



3e, 



(14) 



where we set ey = for tetragonal symmetry. 

We now compare our results with those of Eq. (14) of BS. In that equation the only 
hopping matrix element that was included was that between the ground state orbitals of the 
Cu ions. If we keep only such terms in Eq. (0), we obtain 



rams rams 

XX yy 



— 2A^tni 



35 + C 



00 o 



el{ex + A-2B + Cy 



Px X 



1 



1 



ex + A-5B tx + A + B + 2C 

(15) 



rams 
(7 ~~ 



2 .2 



C 



2+A 



0,Px 



00 9 



e2(e, + A + 45 



g2 g2 



1 



1 



+ A + 45 + A + 45 + 2C, 



(16) 



where we used Eq. ( ^3]) to set too = to,p^/^Px wrote our expression in terms of singlet 
and triplet energy denominators to facilitate comparison with BS. We note the following 
differences between their results and ours. 1) Our results are smaller by an overall factor 
of 2. 2) Instead of our evaluation in which too = tpd/^ (^*^ their notation), they use 
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too = ildl + Exy)- 3) The last energy denominator in their Jyy is wrong: (their 45 should 
be replaced by B). With respect to the first difference we would note that, as described 
in the text, we did compare results from the full diagonalization with those using the spin 
Hamiltonian. Thus an error by a factor of two in our calculations is extremely unlikely. 
Difference 2 comes about because BS do not sum over all processes. In particular consider 
process 1 of Fig. 7. When the hop actually consists of two hops, one from a Cu to an 
O and another from an O to a Cu, this process corresponds to two orderings of the three 
perturbations, one in which the spin-orbit interactions comes first and one in which it comes 
second. (The case when it comes third should be identified with processes 2 or 3 of Fig. 7.) 
Summing over these two orderings converts the denominator {/S. + Exy) of BS into ours. The 
correct energy denominator can also be obtained from Eqs. (El) and (E2). 

In principle, we ought also to compare with Eq. (20) of BS, where contributions involving 
ipi (r) are claimed to be included. Here their results are so different from ours that we can not 
identify their terms with ours. In particular, we note the following. 1) Although they claim 
to include the effects of V'i(i')) their expressions do not include any energy denominators 
which depend on the associated crystal field energy ei. Obviously, when ei = 0, the ground 
manifold would be described by a totally different spin Hamiltonian to remove the spin and 
orbital degeneracy. 2) In our Eq. (0) [ |18| toi enters in several places and thus gives rise 
to many more terms than appear in BS. 3) When toi 7^ 0, as we have noted, each bond 
has biaxial anisotropy in contrast to the axial anisotropy they implicitly assume. 4) We 
include hopping between excited crystal field states. These give contributions which are of 
the same order of magnitude as those involving hopping into the ground state. 5) We have 
not included covalency corrections. It may indeed be a good idea to include such corrections, 
but at present the parameters are themselves so uncertain that we regard this correction as 
a refinement. 
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TABLES 

TABLE I. Values in eV of Parameters Used 



Tja. 












ex = el 






^Py 


4 


9.34 


7.00 


0.15 


0.58 


1.5 


1.8 


1.8 


1.8 


3.25 


3.25 


3.25 



8.8. Local density calculations of Refs. 28 and 2( give 10.5 and 



^ Refs. and g| use [/q 
9.4, respectively. 

^ The Racah coefficients, A, B, and C are defined in Ref. 
For the solid, the values of B and C are appropriately taken from the free-ion optical 



values of Ref. as is discussed by Eskes et al.. [ ^] 
See Ref. M 



' See Ref. gl. 
^ Ref. |32| gives ep 



3.5, but smaller values of ep are plausible. [ |T7| 
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FIGURES 

FIG. 1. A CuO plaquette. Here we distinguish between "y" oxygen ions (on y-directed bonds) 
and "x" oxygen ions (on x-directed bonds.) In each case we show a p orbital on the oxygen ion to 
which we give the symmetry label, z, since these orbitals can only hop to an orbital on a copper 
ion with that same symmetry label, i.e., to tpz ~ xy, which is also shown. 

FIG. 2. Schematic view of dxz and dyz of two Cu ions when they are on the x~axis. Note that 
while dxz orbitals are in the same plane, the plane of dyz orbitals are parallel to each other. Hence 
there is no reason that tx^x should be equal to ty^y. 

FIG. 3. Comparison of perturbation results (dotted line) with the exact results (solid line). 
Here Jy, J±, and Jz correspond to Jxx, Jyy, and Jzz of Eq. (|50|) . The hopping matrix elements 
ta^l3 are estimated from Eq. (|5^ ) as explained in the text. In the left and right panel the ta,i3 and 
Ka,/3 are replaced by t ta,/3 and K Ka,/3 respectively. The values of A (in eV), t, and K are given 
in the panels. Preprint version: labels for curves didn't print well. In each panel the upper curves 
are J± — Jz and the lower curves are J±-J\\ ■ 

FIG. 4. Noninteracting spin-wave spectrum along different high symmetry lines in the Brillouin 
zone according to Eq. (^) for 9 = 0. For this plot the lattice constant a is set equal to unity. The 
J's and uj are all in the same arbitrary units. 

FIG. 5. Variation of the zero-point energy as a function of the angle 9 between the staggered 
magnetization and a [1,0] direction in the easy plane. Dotted and solid lines are from the approx- 
imate [Eq. (^)] and exact [Eq. (pSf)] expressions for Ez{9), respectively. The J's and Ez{9) are 
all in the same arbitrary units. 

FIG. 6. Anisotropies (energy differences between triplet states) as a function of t and A for 
non-constant (left) and constant (right) on site Coulomb interaction. If the energies of the triplet 
states are Ai < A2 < A3, the data points are Ai — A2 (circles) and Ai — A3 (diamonds). The solid 
and dashed lines are power-law fits, as indicated. 

FIG. 7. Processes (1-5) which contribute to Qa^ [in the order written in Eq. (H9)] and those 
(6-7) which contribute to Qa^ [in the order written in (H13)]. Here the left site is the site i and 
the right site is site j. The dashed line depicts the first matrix element and the full line the second. 
The orbitals are the ground state lowest, the state |1) next, and the state \a) highest. In term 3 the 
second process promotes the left-hand hole to an excited state, whereas in term 2 the left-hand 
hole remains in the ground state. 
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